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16 Abstract 

This report investigates several applications of low lift to drag ratio aerobrakes which use angle- 
of-attack variation for control. These applications are: return from geosynchronous or lunar orbit to low 
Earth orbit and planetary aerocapture at Earth and Mars. A number of aerobrake design considerations are 
reviewed. It was found that the flow impingement behind the aerobrake and the aerodynamic heating 
loads are the primary factors that control the sizing of an aerobrake. The heating loads and other loads, 
such as maximum acceleration, are determined by the vehicle ballistic coefficient, the atmosphere entry 
conditions, and the trajectory design. Several formulations for defining an optimum trajectory are 
reviewed, and the various performance indices that can be used are evaluated. The “nearly 
grazing” optimal trajectory was found to provide the best compromise between the often conflicting goals 
of minimizing the vehicle propulsive requirements and minimizing vehicle loads. The relationship 
between vehicle and trajectory design is investigated further using the results of numerical simulations of 
trajectories for each aerobrake application. The data show the sensitivity of the trajectories to several 
vehicle parameters and atmospheric density variations. The results of the trajectory analysis show that 
low lift to drag ratio aerobrakes, which use angle-of-attack variation for control, can potentially be used 
for a wide range of aerobrake applications. 
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TECHNICAL MEMORANDUM 


APPLICATIONS OF LOW LIFT TO DRAG RATIO AEROBRAKES USING 
ANGLE OF ATTACK VARIATION FOR CONTROL 

I. INTRODUCTION 


As the utilization of outer space increases for a wide range of endeavors, from commercial 
ventures to human exploration of the solar system, the search for technologies and techniques which 
would maximize the benefits is increasing. One technological area that is receiving considerable 
attention is space transportation. If space transportation can be performed with maximum efficiency, 
the cost can be minimized. One of the primary objectives in studies of space transportation systems 
is to optimize vehicle performance so that payload mass or operational capabilities can be 
maximized. 

Spacecraft aerobraking is one transportation technology that has great potential for increasing 
the performance capabilities of future space vehicles. Aerobraking is a technique in which a space- 
craft passes through a region of a planet's atmosphere in order to alter its orbital trajectory. During 
the flight through the atmosphere, aerodynamic forces are used to change the spacecraft's velocity 
and direction of flight. The techniques and design issues associated with aerobraking and the 
potential applications of aerobraking have been investigated for many years. It is believed that 
aerobraking offers significant potential performance increases in applications related to orbital 
transfers and planetary orbital capture. 


A. Applications of Aerobraking 

During the early 1960's, the application of aerobraking to orbital transfers was first proposed. 
In 1962, H.S. London established that aerobraking could provide substantial performance advant- 
ages for orbital plane change maneuvers [1]. This application of aerobraking was called the syner- 
getic plane change. This type of maneuver uses aerodynamic forces directed out of the orbital plane 
to turn the vehicle to a new direction. The synergetic plane change was studied extensively during 
the 1960’s with peak activity around 1967. These studies focused on slender reentry vehicles with 
high lift to drag ratios (L/D) and metallic thermal protection systems (TPS) [2]. These concepts 
were the forerunners of the space shuttle and national aerospace plane. 

The synergetic plane change maneuver is illustrated in figures 1-1 and 1-2 [3], The syner- 
getic plane change is accomplished with a combination of propulsive maneuvers plus the aerobraking 
maneuver [2], The first maneuver is a propulsive deorbit to a shallow atmospheric entry. Some of 
the plane change could be performed during this maneuver if it would optimize the overall vehicle 
performance. Once in the atmosphere, the vehicle banks in order to rotate the lifting force out of the 
original orbit plane. This aerodynamic turn results in a rotation of the orbit plane [1]. This is illus- 
trated in figure 1-2 [4], In addition to a change in the orbit plane, the spacecraft's velocity is 
decreased during the aeropass. This results in a lowering of the orbit apoapsis. 

In many applications of aerobraking it is necessary to perform a specified plane change and 
lower the orbit apoapsis to a specified altitude. The plane change that can be accomplished without 
causing the exit apoapsis to fall below the target altitude is often called the "free" plane change. 



DEBOOST AV 


INITIAL ORBIT 


VEHICLE TRAJECTORY 



Figure 1-1. Aeroassisted orbital transfer. 
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Figure 1-2. Atmospheric plane change. 


since no propellant is needed [4]. If the desired plane change is greater than the achievable free 
plane change then a propulsive manuever could be used to perform part of the plane change or raise 
the orbit apoapsis back up to the target altitude. As shown in figure 1-1, this maneuver could take 
place either during the aeropass or immediately after atmospheric exit [1]. 

The final propulsive maneuver for the synergetic plane change occurs when the vehicle 
reaches the target apoapsis. This maneuver raises the orbit periapsis out of the atmosphere to the 
desired final periapsis altitude. In the synergetic plane change, propulsion is used primarily to lower 
and raise the orbit apsides, and aerodynamic force is used primarily to change the orbit plane [1]. 

Another application of aerobraking to orbital transfers makes use of aerodynamic forces for 
velocity changes only. In this application there is less need for lift since there is no atmospheric 
plane change. This type of aerobraking maneuver is called a coplanar aeroassisted orbital transfer. 
The purpose of the atmospheric pass is solely to dissipate energy in order to lower the apoapsis 
altitude [1]. This aerobrake application was studied in the 1960’s and again in the early 1980’s. 

These studies led to innovative approaches for carrying cargo to geosynchronous orbit or the Moon 
from a space station in low Earth orbit (LEO) using vehicle concepts with low L/D and low ballistic 
coefficients. They indicated that the application of aerobraking to orbital transfers could increase 
vehicle payload capabilities by as much as a factor of 2 over conventional all propulsive transfers [2]. 
The mission profile for a coplanar aeroassisted orbital transfer is similar to the synergetic plane 
change except there is no plane change in the atmosphere and there is no apoapse raising maneuver 
in the atmosphere or at atmospheric exit (fig. 1-1). 

The significant potential for vehicle performance gains for aeroassisted orbital transfers in 
Earth orbital applications also exists for interplanetary missions. The application of aerobraking to 
planetary orbit capture maneuvers is an attractive option which was first studied from 1964 to 1968 
and again beginning in 1979 [2], 

Planetary aerocapture is illustrated in figure 1-3. A spacecraft arrives at a planet on a 
hyperbolic trajectory (relative to the planet) which penetrates the atmosphere. Once the vehicle 
enters the atmosphere, the aerodynamic drag slows the vehicle to a desired orbital velocity. Ideally, 
the velocity is reduced so as to achieve a target parking orbit apoapsis. When the vehicle reaches 
the apoapsis, a propulsive maneuver is used to raise periapsis out of the atmosphere to its final 
value. Just as in the coplanar aeroassisted orbital transfer, the purpose of the atmospheric pass is 
primarily to dissipate energy, however, there is no reason that a plane change could not also be 
performed. Aerocapture is essentially a transfer from a hyperbolic orbit to an elliptical orbit. 


B. Topics of Discussion 

There are numerous types of problems that must be addressed in the study of aerobrake 
applications. These problems include thermodynamics of hypersonic flight, TPS options, L/D selec- 
tion and operational considerations such as control methods and guidance strategies. The solution of 
these problems is not always straight-forward and usually requires numerical methods. This reliance 
on numerical solutions which depend on many parameters makes it difficult to obtain truly general 
results for aerobraking analyses. The theory of optimal aeroassisted orbital transfer is considered to 
be in only the development stage. The goal of this development is to establish a sound basis for 
determining the best strategies for utilizing aerobraking and to describe the characteristics of the 
optimal strategies [1]. This report is an attempt to characterize the techniques and performance 
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of one aerobraking strategy as it could be applied to a number of applications. The strategy that will 
be studied will be the use of low L/D aerobrakes with angle of attack control. To better understand 
the reason for the selection of this strategy it is helpful to consider some of the topics which are 
being addressed in the current studies of aerobraking. 

Two interrelated topics which are being addressed in aerobrake studies are aerobrake L/D 
requirements and aerobrake flight control strategies. Early studies of aerobraking focused on high 
L/D concepts while more recent studies have considered low L/D concepts. The high L/D aerobrakes 
have definite advantages in performing plane changes. However, since the lift force does not aid in 
slowing the vehicle down, there is no real need for high L/D for many orbital transfer applications. 
Low L/D aerobrakes may also offer more efficient payload packaging and better protection from the 
aero thermal environment than high L/D aerobrakes. If there is no atmospheric plane change, lift force 
is needed only for vehicle control. The amount of control required is dependent on the atmospheric 
conditions encountered, vehicle ballistic coefficient, and the control strategy that is employed. In this 
report the performance characteristics of low L/D aerobrakes for several mission applications will be 
determined, and some of the requirements for adequate vehicle control will be investigated. 

The aerobrake control strategy is closely related to the L/D selection. Almost all aerobrake 
concepts that have been investigated have assumed "roll control" to guide the vehicle through the 
atmosphere. In this control strategy, ascent and descent in the atmosphere is controlled by rotating 
the lift vector about the velocity vector (bank angle). The bank angle is varied in order to change the 
vertical component of the lift vector and achieve a desired rate of ascent or descent. Any unwanted 
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lateral component of lift is cancelled by alternating the bank angle from side to side. One reason that 
roll control has been assumed in many studies is because it was successfully used for reentry control 
for the Mercury, Gemini and Apollo spacecraft. It will also be used for the Aeroassisted Flight 
Experiment. An alternative control strategy is to vary the magnitude of the lift vector by rotating the 
vehicle about the pitch axis. This strategy varies lift by changing the vehicle angle of attack. It has 
not been studied as much as roll control but offers several potential advantages such as faster 
response times and elimination of lateral forces. 

The second section of this report will address aerobrake design considerations. First, the 
aerodynamic characteristics of aerobrakes will be discussed. This discussion will begin with a 
description of the aerodynamic coefficients, followed by an explanation of the conditions required for 
static stability. The significance of L/D and ballistic coefficient will also be described. Another aspect 
of aerobrake design that will be discussed is aerodynamic heating. A discussion of aerobrake con- 
struction will address thermal protection system design and aerobrake mass estimation. 

Section III will discuss the optimization of aerobrake trajectories. This will include two 
formulations that can be used, depending on the choice of a performance index. Several potential 
performance indices will be described and evaluated. The process used to determine a lift control 
function will be outlined. This control function will be used as a basis for the angle-of-attack vari- 
ation for the aerobrake applications investigated in this report. To complete the discussion of the 
optimum control function, several vehicle control strategies which could be used will be described. 

Section IV will contain the bulk of the data that was generated. This chapter will begin with a 
description of the computer program that was used to simulate the aerobrake trajectories and the 
techniques used to model the atmosphere. Data will include constant and variable angle-of-attack 
trajectories for several aerobrake applications. These applications include return from GEO, return 
from Lunar orbit. Earth aerocapture, and Mars aerocapture. The data will show the sensitivity of 
aerobrake trajectories and vehicle loads to angle of attack, ballistic coefficient, and atmospheric entry 
conditions. It will also be used to define target periapsis tolerances. 

The final section of this report will summarize the characteristics of low L/D aerobrakes and 
the results of the trajectory and performance analyses. This section will include an overall assess- 
ment of the low L/D aerobraking strategy and summarize the conclusions that have been reached 
concerning aerobrake control requirements. 

The study of aerobraking and its application to future space missions involves the analysis of 
a number of complex interrelated issues. As mentioned earlier, there is no definitive theory that can 
be applied to the analysis of aerobraking. Instead, current studies seek to identify the potential 
strategies for the application of aerobraking and the methodology to evaluate these strategies. The 
objective of this report is to study one possible control strategy for aerobraking and to evaluate the 
potential benefits for a number of applications. It is hoped that the results obtained and the analysis 
performed in this report will add to the overall body of knowledge from which a general theory of 
aerobraking optimization may be derived. 


II. AEROBRAKE DESIGN CONSIDERATIONS 


This section will describe some of the important design considerations that must be made in 
the design of aerobrakes. The most important design considerations for aerobrakes arise from their 
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aerodynamic characteristics and the flight environments in which they operate. The unique con- 
figurations of aerobraked vehicles give them aerodynamic characteristics that are quite different from 
other atmospheric vehicles. These aerodynamic characteristics influence the sizing of the aerobrake 
and the stability of the vehicle during flight. The aerodynamic properties also determine the trajec- 
tory requirements for a particular orbital transfer and the dynamic and thermal loads experienced by 
the vehicle. These loads are also important factors which determine the mass of the aerobrake and 
the performance gain that is possible. 


A. Aerodynamic Characteristics of Aerobrakes 

1. Aerobrake Lift and Drag 

There is a degree of uncertainty that exists in the analysis of the aerodynamic characteristics 
of aerobrakes since aerobraking is performed in the upper atmosphere. Figure 2-1 shows an 
approximate altitude profile of an aerobraking trajectory and the flow regimes that are experienced 
[5]. An aerobraked vehicle spends considerable time in the noncontinuum flow regimes [5]. Unfor- 
tunately, there has been little research in the aerodynamics of thermodynamics of noncontinuum 
flows [6]. Most of the aerodynamic data for aerobrakes has been developed using Newtonian aero- 
dynamics for hypersonic conditions. Under these conditions, the dependence of the aerodynamic 
coefficients on Mach number and Reynolds number is disregarded [3], therefore, the aerodynamic 
coefficients depend only on angle of attack. 

Many potential aerobrake shapes have been investigated. In this report, only low L/D aero- 
brakes will be described. Figure 2-2 shows three typical low L/D aerobrakes. The spherical aero- 
brake is based on the heat shield of the Apollo command module. The use of this shape allows the 
application of past research and flight data to aerobraking analyses. The sphere-cone configuration 
consists of a circular cone with a spherical section at the tip. The raked-off elliptical cone is defined 
as a cone with an elliptical cross section and a spherical tip which is "raked-off at an oblique angle 
as shown in figure 2-3 [7]. This shape provides a higher L/D at low angles of attack. The con- 
figuration shown will be used for the planned aeroassist flight experiment (AFE). 

Figure 2-4 shows the lift and drag coefficients for the three aerobrake shapes [8]. The aero- 
dynamic coefficients for the spherical (Apollo), sphere-cone, and AFE configurations were calcu- 
lated using the Mark IV Supersonic-Hypersonic Arbitrary Body Program [9]. This program uses the 
modified Newtonian pressure method to calculate pressure coefficients for each element of surface 
area on the aerobrake. The components of the pressure coefficients are then used to calculate the 
aerodynamic coefficients. The wind tunnel data shown in figure 2-4 was obtained using a 4-in 
diameter aerobrake model at a Mach number of 10.2 [10], The AFE has been designed to fly at a 
specific angle of attack of 17° which puts the stagnation point at the center of the spherical tip. The 
aerodynamic data at other angles of attack are shown for illustrative purposes only. 

The curves shown in figure 2-4 give an indication of the unique aerodynamic characteristics of 
aerobrakes. The drag decreases as the magnitude of the angle of attack increases. This is opposite 
of most aerodynamic vehicles. Another characteristic that is opposite of most aerodynamic vehicles 
is that aerobrakes experience positive lift (lift up) at negative angles of attack and negative lift (lift 
down) at positive angles of attack. This characteristic can be used to differentiate between aero- 
brakes and aerobodies (of lifting bodies). Vehicles having positive lift at positive angles of attack 
can be defined as aerobodies [11]. For this definition to hold true for nonsymmetrical aerobrakes 
such as the AFE, a further stipulation must be made that the bank angle is zero. If the AFE is flown 
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inverted (bank angle = 180°) then the lift is positive at positive angles of attack. The advantage of 
this will be explained later. It should be noted that there is no standard definition of aerobrake 
geometry or coordinate systems so care must be taken when referring to angles of attack and bank 
angles. 


The direction of lift is an important factor in aerobraking because it influences the vertical 
movement of the vehicle in the atmosphere. Lift should work in harmony with drag in controlling the 
trajectory of the vehicle. If more drag is needed, the lift force should pull the vehicle down (deeper) in 
the atmosphere. If less drag is needed, the lift force should move the vehicle upward in the atmos- 
phere [12]. Referring back to figure 2-4, it can be seen that this coordination of lift and drag exists at 
negative angles of attack for symmetrical aerobrakes (spherical, sphere-cone) and at a bank angle of 
180° for the AFE (raked-off cone). Even though the most efficient coordination of lift and drag is 
achieved when lift is directed upward, it will be shown later that when a trajectory is optimized to 
minimize certain vehicle loads it is sometimes necessary to fly the aerobrake with lift directed 
downward. 


2. Aerobrake Static Stability 

One of the most important design considerations for aerobraked vehicles is the aerodynamic 
stability of the vehicle. The vehicle must be designed to maintain the proper flight attitudes in order 
to achieve the desired trajectory. It is beyond the scope of this report to investigate the dynamic 
stability of aerobraked vehicles, therefore, the discussion of aerodynamic stability will be limited to 
the necessary conditions for vehicle equilibrium (trim) and longitudinal static stability. These are 
preconditions for dynamic stability. 

In the discussion of vehicle trim, it is necessary to define a vehicle coordinate system. This 
coordinate system is shown in figure 2-5. For the present discussion only, the spherical (Apollo) 
aerobrake shape will be considered. The x axis is defined by the central axis of the vehicle. The z 
axis is defined by the longitudinal plane of the vehicle. The origin is defined to be the center of the 
front of the aerobrake. In defining the aerodynamic forces and moments, the reference area (A) is the 
cross-sectional area of the aerobrake, the reference length (d) is the aerobrake diameter, and the 
moment reference point is the origin (jc = 0, z = 0). It should be noted that the location of the moment 
reference point is arbitary, but it is tied to a particular set of aerodynamic coefficients. The aero- 
dynamic coefficients for a spherical aerobrake are listed in table 2-1 [8], 

When the vehicle attitude is maintained in a particular orientation, the aerodynamic forces 
and moments are in equilibrium, and the vehicle is in a "trimmed" configuration. The necessary 
condition for trim is that the sum of the moments about the vehicle center of gravity (eg) is zero. If 
the vehicle eg is located at the arbitary position shown in figure 2-5, the moments can be summed as 
follows: 


M C g — Cm cj clA + C,y (X C g—XR)q^-CA(Zcg—ZR)q^ ; (2-1) 

the moment coefficient about the eg can be defined as: 


_ A/ Cg 

'-"leg - , . 

q d A 


( 2 - 2 ) 
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Figure 2-5. Aerobrake coordinate system. 


Table 2-1. Aerodynamic coefficients for a spherical aerobrake. 

REFERENCE AREA = AEROBRAKE FRONTAL AREA 

REFERENCE LENGTH = AEROBRAKE DIAMETER 



MOMENT REFERENCE 

POINT = 

FRONT CENTER 

OF AEROBRAKE 

X 

ii 

■< 

it 

N 

II 

o 

b 


ALPHA 

L/D 

CL 

CD 

CA 

CN 

CM 

(DEG.) 







0 

0.0 

0.0 

1.5167 

1.5167 

0.0 

0.0 

4 

-0.0602 

-0.0908 

1.5072 

1.5098 

0.0146 

-0.0135 

6 

-0.0905 

-0.1353 

1.4953 

1.5013 

0.0218 

-0.0202 

8 

-0.1 208 

-0.1787 

1.4789 

1.4894 

0.0289 

-0.0268 

1 0 

-0.1514 

-0.2207 

1.4580 

1.4741 

0.0359 

-0.0333 

1 2 

-0.1821 

-0.2609 

1.4328 

1.4557 

0.0427 

-0.0396 

1 4 

-0.2131 

-0.2990 

1.4034 

1.4340 

0.0494 

-0.0457 

1 6 

-0.2443 

-0.3347 

1.3702 

1.4094 

0.0559 

-0.0516 

1 8 

-0.2759 

-0.3678 

1.3334 

1.3818 

0.0622 

-0.0572 

2 0 

-0.3078 

-0.3980 

1.2933 

1.3514 

0.0683 

-0.0626 

22 

-0.3400 

-0.4251 

1.2501 

1.3183 

0.0742 

-0.0677 

24 

-0.3727 

-0.4489 

1.2043 

1.2828 

0.0798 

-0.0724 

26 

-0.4058 

-0.4693 

1.1562 

1.2449 

0.0851 

-0.0768 

28 

-0.4394 

-0.4860 

1.1062 

1.2049 

0.0902 

-0.0809 

3 0 

-0.4733 

-0.4992 

1.0546 

1.1629 

0.0950 

-0.0846 
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Using this definition, equation (2-1) can be rewritten in the following form: 


n — n , t" <*c g Xr) n (Z cg Z/j) 

Cm “ - c m +c^ - c A - • {2 . 3 ) 

The condition for trim is: 

M C g = C mc g = 0 . (2-4) 

In deriving equation (2-3), an arbitrary eg location was used. The combinations of ;c cg and z C g 
that satisfy equation (2-3) define a line in the longitudinal plane called the trim line. This line can be 
defined by setting equation (2-3) equal to zero then solving for z cg (normalized to z C gfd) 

Zcg _ C m Qf Xjt | Zft | Cat ^cg 

d Q Ca d d D a d (2-5) 


z intercept slope 

As long as the vehicle eg falls on the trim line, the vehicle will maintain a constant attitude. Since C m , 
Cn, and Ca are functions of angle of attack, the trim line is also a function of angle of attack. Figure 
2-6 shows the trim lines at various angles of attack for a spherical aerobrake [13]. 

Placing the eg on the desired trim line insures that the vehicle will be in equilibrium, but it 
does not guarantee that it will be in a stable equilibrium. Longitudinal static stability depends on the 
change in pitching moment when a disturbance causes a slight change in the angle of attack. The 
necessary condition for longitudinal static stability is defined as: 


dC/rkt 

da 


= C 


n ta 


<0. 


( 2 - 6 ) 



Figure 2-6. Trim lines for a spherical aerobrake. 
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Figure 2-7 shows an example of an aerobrake trimmed at an angle of attack of 15°. The plot 
of pitching moment coefficient versus angle of attack shows that C OT( a is negative and C mcg = 0 at a = 

15°. These conditions indicate that the vehicle possesses longitudinal static stability. For any 
disturbance causing an increase in angle of attack (clockwise rotation), the aerodynamic forces will 
cause an increase in Cm C g in the negative direction (counter-clockwise rotation). The opposite is true 

if a disturbance causes a decrease in angle of attack [12]. 



CONDITION FOR STATIC STABILITY: (W<0 

FOR ANY DISTURBANCE CAUSING AN INCREASE IN cc (CLOCKWISE ROTATION) 

THE AERODYNAMIC FORCES WILL RESULT IN AN INCREASE IN C m IN THE NEGATIVE DIRECTION 
(COUNTER-CLOCKWISE ROTATION) 

Figure 2-7. Longitudinal static stability during aerobraking. 

The value of C ma depends on the aerodynamic coefficients and the eg location. The relation- 
ship can be found by differentiating equation (2-3) with respect to angle of attack then using equa- 
tion (2-5) to eliminate z c g from the result. 

= dC m bC A Cjn , (dCv dCt Qy\ Ucg -*/?) 
m “ da da C A Ucc da C A I d ' (2-7) 


For a given trim condition (angle of attack), the aerodynamic coefficients are constant so the stability 
depends on the eg location. The eg location at which C ma = 0 is called the neutral point. The x coor- 
dinate of this point can be found by setting equation (2-8) equal to zero and solving for x cg : 


c A 


dCm 

dC A 


&np 


dCfj_ 

dC A 


C A 


+ X R . 


( 2 - 8 ) 
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For the example in figure 2-7, at a trim angle of attack of 15°, the neutral point is at x = 0.7284 d. If 
the jr C g is greater than x np , then C OT(X will be positive and the vehicle will be unstable. Obviously it is 
desirable to have the eg as far in front of the neutral point as possible. Equations (2-7) and (2-8) 
show that the stability is also a function of aerobrake diameter. As the diameter increases, C m<x 
becomes negative, and the x np increases so the vehicle becomes more stable. This indicates that 
stability considerations should be included in the process of sizing an aerobrake. 


3. Flow Impingement Behind Aerohrakes 

One of the more basic considerations in sizing an aerobrake is the flow condition behind the 
aerobrake. As the flow passes the edge of the aerobrake, it turns inward toward the center of the 
vehicle. The aerobrake must be sized such that this flow does not impinge on the payload behind the 
aerobrake. Figure 2-8 shows the impingement angle behind a spherical aerobrake as a function of 
angle of attack [14]. The data was extracted from reference 20. The most critical impingement occurs 
at the most forward edge of the aerobrake. The flow impingement angle is equal to the sum of the 
magnitude of the angle of attack and the flow turning angle. 



ANGLE OF ATTACK, a (DEG) 

Figure 2-8. Flow impingement angle for an aerobrake. 
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4. Lift to Drag Ratio 

It was indicated that the combined effects of lift and drag are important in controlling an aero- 
brake trajectory. This relationship has led to the use of the L/D as an important parameter in charac- 
terizing aerobrakes. The LTD, along with the ballistic coefficient, determines the aerodynamic control 
capability and influences the heating and acceleration loads on the vehicle [1]. Figure 2-9 shows the 
L/D and ballistic coefficients for typical aerospace vehicles [2]. L/D’s greater than 1.5 require 
blended wing-body configurations. Most aerobrakes fall into the low- to mid-L/D range. The primary 
factor in choosing an L/D for an aerobrake is the type of mission to be flown. High-L/D aerobrakes 
(L/D > 2.0) have significandy greater plane change capability in the atmosphere than low L/D aero- 
brakes (L/D < 0.8). Unfortunately, the benefit of using a high L/D aerobrake to perform a synergetic 
plane change is counteracted by an increase in structural mass. In comparing high L/D aerobrakes to 
low L/D aerobrakes for missions returning from GEO, Talay et al. [2] found very little difference in 
payload capability, but the high L/D aerobrake had significantly higher aerodynamic heating and 
structural mass. These results were verified by Mease and Vinh [15]. Table 2-2 shows the relation- 
ship between L/D and several trajectory parameters for an aerobraked transfer from GEO to LEO. 
Based on these results, it appears that high L/D offers little or no advantage for this mission [2], 

The previous conclusion could be interpreted from a different perspective as pointing out a 
disadvantage of the synergetic plane change rather than high L/D. The reason that the thermal and 
structural requirements are greater for the high L/D configuration is that the vehicle trajectory went 
deeper into the atmosphere to maximize the synergetic plane change. In aerobrake applications with 
no synergetic plane change, it has been found that increasing the L/D results in a decreased penetra- 
tion into the atmosphere (if lift is downward) which decreases the vehicle loads [14). Unfortunately 
in these cases, the mass penalty associated with the complex aerodynamic shapes can make high 
L/D configurations unattractive. It appears that high L/D aerobrakes do not provide any significant 
advantages for missions to high Earth orbit such as GEO or for missions that do not require syner- 
getic plane changes such as return from lunar orbit or planetary orbital capture. This is one of the 
primary reasons for considering only low L/D aerobrakes in this report. 

Figure 2-10 shows the L/D as a function of angle of attack for low L/D aerobrakes [8]. The 
symmetrical aerobrakes have very similar L/D characteristics while the raked-off cone (AFE) L/D is 
slightly higher in magnitude. Based on the analysis performed by Mayo, Lamb, and Romere [7], the 
L/D for low L/D aerobrakes at low angles of attack can be approximated by the following relation- 
ship: 


L. 

D 


= sin (a) . 


(2-9) 


This formula is useful in making a first guess for the angle of attack required to achieve a desired 
L/D. 


The choice of an L/D for an aerobrake application is complicated by the uncertainty in the 
aerodynamics of the vehicle in the upper atmosphere. G.D. Walberg has indicated that the L/D 
decreases at high altitudes due to viscous effects [2]. This phenomenon has been ignored in most 
aerobrake studies [1]. Maslin was the first to analyze aerobrake trajectories with this effect 
included. He developed an approximation for the L/D degradation at high altitudes [2]. Figure 2-11 
shows the variation of maximum L/D with altitude for the space shuttle [2]. The solid curve 
represents wind tunnel and flight data, and the hatched band represents Maslan's approximation. 
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Figure 2-9. L/D's and ballistic coefficient values for typical 
lifting bodies [2], 


Table 2-2. Aerobrake parameters for a GEO to LEO transfer. 

LIFT TO DRAG RATIO 



LOW 

MODERATE 

HIGH 

MAXIMUM L/D 

0.845 

1.5 

2.9 

MINIMUM ALTITUDE (KM) 

58.8 

58.2 

51.5 

MAX. DYNAMIC PRESSURE 
(N/M**2) X 0.001 

1 5.9 

18.6 

44.2 

MAX. CONVECTIVE HEATING 
RATE (W/CM**2) 

193.1 

222.8 

361.4 

MAX. ACCELERATION 
(EARTH g'S) 

3.6 

2.7 

1.8 
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Figure 2-10. L/D's for aerobrakes. 





While it appears that Maslin's approximations underpredict the L/D at high altitudes, the shuttle 
data shows that there is still significant degradation of the L/D above 70-km altitude. 


5. Ballistic Coefficient 

The ballistic coefficient is a very important design parameter for aerobrakes. As mentioned 
earlier, it is closely related to the L/D in its influence on the vehicle and trajectory. The ballistic 
coefficient determines the structural and heating loads that are experienced during the aeropass as 
well as several other trajectory parameters. One of the objectives of this report is to investigate the 
relationship between the ballistic coefficient and the aerobrake control requirements. This relation- 
ship will be described in subsequent sections. 

The ballistic coefficient is defined by the following equation: 

p = -2t_. 

C d A ( 2 - 10 ) 

It provides a measure of the capability of an aerodynamic drag force to change the kinetic energy 
(velocity) of the vehicle. The mass is an indicator of the vehicle's inertia, and the product of the drag 
coefficient and aerobrake cross-sectional area is an indicator of the drag force that can be generated. 
A high ballistic coefficient indicates that a larger drag force is required to achieve a given velocity 
change. For a given aerobrake ballistic coefficient, an increase in drag is achieved by flying the 
vehicle deeper in the atmosphere to achieve a greater dynamic pressure. This increases the vehicle 
structural and thermal loads [14]. 

Figure 2-9 shows that for most aerobraked vehicles the ballistic coefficients will not exceed 
500 kg/mA Figure 2-12 shows the sensitivity of the ballistic coefficient to vehicle mass and 
aerobrake diameter. It shows that there are extremes to be avoided. There is a point at which 
increasing the diameter of the aerobrake has little benefit in lowering the ballistic coefficient due to 
the increase in aerobrake mass. There is also a point at which decreasing the diameter of the 
aerobrake (to reduce the mass) will result in a greatly increased ballistic coefficient which could 
result in an increase in structural and thermal loads.. 


6. Aerodynamic Heating 

Aerodynamic heating is probably the most critical design factor that must be considered in 
aerobrake design. The thermal loads that are anticipated in most aerobrake applications require 
insulation materials which are at the limit of current technology. The heating loads influence the 
aerobrake material selection and sizing (diameter) and also constrain the atmospheric trajectories 
that can be flown. Unfortunately, the understanding of the thermal environments that will be 
encountered by aerobraked vehicles is not complete. There is evidence that chemical and thermo- 
dynamic nonequilibrium caused by the dissociation and recombination of air molecules will signifi- 
cantly influence the convective and radiative aerodynamic heating [2], One of the primary objectives 
of the planned AFE is to determine the extent of nonequilibrium effects. There are two principle 
phenomena that may be present. They each have a different effect on the heating conditions. 

Chemical nonequilibrium is the result of the catalysis of the aerobrake surface. The high 
temperatures near the aerobrake cause the air molecules to dissociate. A certain fraction of the 
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atoms that strike the aerobrake surface would recombine due to catalysis effects. When they 
recombine, the energy of dissociation is left on the wall of the aerobrake and contributes to the heat 
flux [5]. There is evidence that finite wall catalysis could reduce the convective heating rate by up to 
one half of their equilibrium values [2]. Figure 2-13 [5] shows the effect of wall catalysis on the 
peak temperature for an aerobraked return from GEO. 


ANGLE OF ATTACK = 15 DEGREES 



Figure 2-12. Sensitivity of ballistic coefficient to vehicle mass and aerobrake diameter. 



Figure 2-13. Representative peak temperature dependence on catalysis for 
GEO to LEO aerobraked transfers [5]. 
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Another phenomenon that may exist in the aerothermal environment of an aerobrake is non- 
equilibrium shock layer radiation [2]. There is considerable disagreement on the extent of radiative 
heating and its effect. If nonequilibrium radiation occurs, it could increase the peak heating rate by 10 
percent [5] to 100 percent [2], The existence of nonequilibrium radiative heating may be the single 
most important question that exists concerning aeroassisted orbital transfers [2]. This question can 
only be answered through an actual flight test [2]. 

The potential effects of nonequilibrium conditions on the aerodynamic heating are important 
considerations in aerobrake design. Since aerobrakes will have velocities that are greater than 
reentry vehicles, catalysis effects will be even more evident. In the Earth's atmosphere there will be 
substantial nitrogen dissociation, so it would be beneficial to choose insulation materials with low 
nitrogen recombination rates [16], Mars aerobrakes may require different materials to minimize the 
recombination of Mars atmospheric constituents (mostly carbon dioxide). The effects of catalysis 
and convective heating are reduced at higher altitudes due to the decrease in density. This indicates 
that it is desirable to use high-drag aerobrake configurations at very high altitudes [16]. 

In performing aerobrake trajectory analysis, it is impossible to accurately predict the aero- 
dynamic heating environment without the use of Navier-Stokes solutions, viscous shock-layer solu- 
tions, and experimental verification [5]. To avoid these complications, the aerodynamic heating is 
characterized by an indicator that is much easier to calculate. The most common heating indicator is 
the stagnation point convective heating rate, defined by the following equation [1]: 


Q = c p 0 5 V 3 08 W/cm 2 , 


( 2 - 11 ) 


where 


c = 1.83 x 10' 8 r„-°- 5 (1 - gw) 
r n = nose radius (m) 
g w = wall enthalpy/total enthalpy 
p = density (kg/m 3 ) 

V = velocity (m/s). 

The value of the stagnation point heating rate is often used to determine the type and thickness of 
thermal insulation material for an aerobrake, even though wind tunnel tests have shown that the 
maximum heating actually can occur at the edge of the aerobrake [17]. 

For a given aeroassisted orbital transfer, the depth of penetration into the atmosphere 
depends on the L/D and the ballistic coefficient. Since the depth of penetration also influences the 
heating rate, there is a relationship between L/D, ballistic coefficient, and heating rate. This 
relationship for a GEO to LEO transfer is illustrated in figure 2-14 [5]. Note that in this graph the 
heating indicator is the product of the stagnation point heating rate and the nose radius. These data 
indicate that the heating rate is more sensitive to the ballistic coefficient than the L/D [5], 
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Ballistic Coefficient (Kg/m 2 ) 


Figure 2-14. Correlation of peak reference heat flux with ballistic coefficient and 

L/D for GEO to LEO trajectories [5]. 


B. Aerobrake Construction 

1. Aerobrake Thermal Protection 

The aerodynamic heating environments that are anticipated for most aerobraking applications 
will require unique insulation materials to protect the structure of the aerobrake and the payload. 
Currently, there are two general types of thermal protection system materials that may be applicable 
for aerobrakes. 

One type of TPS is the metallic standoff system. This system is composed of a metallic outer 
surface and a layer of insulation between the outer surface and the aerobrake structure. The metallic 
surface radiates some heat away, the rest of the heat is absorbed by the layer of insulation [1]. The 
type of metal used for the outer surface depends on the peak temperature that is experienced during 
the aeropass [18]. Figure 2-15 [18] shows the temperature ranges and weight per unit area for 
several metals. The thickness of the insulation layer depends on the accumulated (or integrated) 
heat load, Q [18], defined by the following equation [1]: 

«.[V 

"0 ( 2 - 12 ) 
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Maximum Surface Temperature 
as a Function of Metallic 
Panel Unit Weight 


Metallic Panel 
Unit Weight lb/ft 2 



Maximum Surface Temperature, T Max - °F 

Figure 2-15. Metallic TPS weight per square foot versus maximum surface temperature [18]. 


To minimize the mass of the metallic standoff thermal protection system, it is desirable to 
have both a low () max and a low Q. Since these quantities have different sensitivities to the trajec- 
tory, a tradeoff between £) max and Q is necessary. Garcia and Fowler found that the total TPS mass 
is reduced by minimizing £) max . The savings gained by reducing the mass of the metallic outer surface 
more than compensates for the increased insulation mass required due to the increased Q [18]. The 
sensitivity of Qm ax and Q to the trajectory, specifically the depth of penetration into the atmosphere, 
will be discussed further in subsequent sections. 


A second type of TPS is called reusable surface insulation (RSI). This type of TPS usually 
consists of ceramic tiles which are coated with a reaction cured glass. The tiles are attached to the 
vehicle structure with a felt strain isolation pad to allow deformation of the skin without damage to 
the tile [14]. This is the type of TPS used on most of the space shuttle. The thickness of RSI for an 
aerobrake depends on the integrated heat load [ 1 8] and the temperature limitations of the material 
[14]. Current RSI materials can withstand a maximum temperature of 1,645 K (2,500 °F) [18]. 

Future materials may have a limit as high as 2,500 K (4,100 °F) [6], which could significantly reduce 
the thermal protection system mass. 

Garcia and Fowler investigated the application of metallic standoff and RSI TPS's for the 
space shuttle [18]. They found that for temperatures less than 1,420 K (2,100 °F) the metallic 
standoff system had a lower mass than RSI. This is indicated in figure 2-16 [18]. Unfortunately, 
most aerobrake applications would exceed this temperature, so reusable surface insulation will be 
required. One possible method for decreasing the thermal load on an aerobrake is to decrease the 
ballistic coefficient of the vehicle by increasing the diameter of the aerobrake [16]. Figure 2-17 [6] 
shows that if ballistic coefficient is reduced to extremely low values, a metallic TPS could be used. 
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Stagnation Heating Rate (w/cm^) MctxifTlUm PsriGl T6mpGr3tUr8 


Metallic TPS 



TPS Weight (Pound x .001 ) 


Figure 2-16. TPS weights for space shuttle entry [18]. 
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Figure 2-17. Sensitivity of thermal loads to ballistic coefficient [6]. 



The benefits of using the lighter TPS would have to be traded against the mass penalties of a larger 
aerobrake. Wurster and Eldred [1] found that a hybrid system that uses both types of TPS may 
provide benefits in reducing the mass and increasing the durability of the aerobrake. 


2. Aerobrake Sizing 

All of the design considerations that have been described influence the sizing of the aero- 
brake. Many tradeoffs are required to achieve the best aerobrake design for a particular mission. The 
ultimate goal in sizing an aerobrake is to minimize the mass. Figure 2-18 [5] shows the aerobrake 
specific mass as a function of accumulated heat load for an aerobrake. The structure mass includes 
the aluminum skin, stiffeners, and frame. The remainder of the mass is the TPS, which consists of the 
tiles, coating, strain isolation pad, and adhesive [5]. Since the mass of the aerobrake depends more 
on the supporting structure than on the insulation, it is advantageous to minimize the diameter of the 
aerobrake. 



In most cases, the minimum diameter is determined by the heating rate constraint of the 
insulation material. The aerobrake must be sized to have a sufficiently low ballistic coefficient such 
that the maximum heating rate will not exceed the limit of the TPS. Figure 2-19 [5] shows the 
relationship between aerobrake diameter, payload mass, and maximum heating rate for a GEO to 
LEO transfer. The specific mass ratio, p b/Cd, corresponds to a value similar to the space shuttle. 

The factor/ represents the deviation of the aerobrake heating rate from a reference sphere. It 
includes the effects of nonequilibrium flow, catalysis, low atmospheric density, and the vehicle shape. 
The value of/ for a spherical aerobrake was derived from data in reference 27. As the payload 
increases, aerobraking becomes more efficient as indicated by a decrease in the aerobrake mass 
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© Raked-Off Elliptical Cone, AFE (t - .9) 
© Spherical Aerobrake, Apollo (f = 1.41) 



Aerobrake Diameter (m) 

Figure 2-19. Vehicle sizing curve based on convective heat flux [5]. 

fraction. Also, if the maximum heating rate for an aerobrake could be increased, the diameter could be 
decreased for a given payload. This, too, would result in a lower mass fraction which would increase 

the efficiency. 

Similar aerobrake sizing data was developed for a lunar transfer vehicle by Dickerson [19]. 
These data were used to construct figure 2-20, which shows several aerobrake parameters as a 
function of the aerobrake payload mass. The aerobrake sizing is based on a maximum heating rate o 
51 W/cm 2 . The values of ballistic coefficient and diameter represent the minimum values that will 
ensure that the heating rate constraint is not exceeded for an aeroassisted return from lunar orbit. In 
some cases, the diameter may need to be increased to prevent flow impingement on the payload. 
Increasing the diameter of the aerobrake decreases the ballistic coefficient and the thermal loads, 
therefore, some of the mass penalty could be reduced by decreasing the insulation thickness. 

The size of the aerobrake appears to depend on two factors, aerodynamic heating and flow 
impingement on the payload. In general, higher L/D and lower ballistic coefficients reduce the vehicle 
loads. Ordinarily, reducing the loads on a vehicle results in a decrease in mass, but this is not 
always the case with aerobrakes. Decreasing the loads may result in a larger than necessyy aero- 
brake, a more complex shape, or worse flow impingement (due to a higher angle of attack) [14], The 
process of sizing an aerobrake and estimating its mass involves a number of trade-offs and Aera- 
tions in which all the design considerations that have been described must be included. Figure 2-21 
[14] shows a schematic of the process used to size an aerobrake. The process can be summarized 
by the criteria listed in table 2-3. In the next section, it will be seen that this aerobrake design 
process must also include careful design of the atmospheric trajectories that will be flown. The 
vehicle design must be closely coordinated with the trajectory design. 
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Ballistic 

Coefficient 


(kg/m 2 ) 

120 - 

118 - 

116- 

114- 

112 - 

110 - 

108 J 


A/B 
Diameter 
(m) 

16 - 

15- 

14- 

13- 

12 - 


A/B 

Mass 

Fraction 

.180- 

. 175 - 

. 170 - 

. 165 - 

. 160-1 



Figure 2-20. Lunar return aerobrake sizing. 



Aerobrake is sized either 
by Impingement or Q 


Figure 2-21. Aerobrake sizing approach schematic [14]. 
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Table 2-3. Aerobrake sizing criteria. 


- GN&C determines L/D requirement 

- L/D determines angle of attack requirement 

- Angle of attack determines impingement angle 

- Ballistic coefficient determines: 

Maximum heating rate 
Maximum dynamic pressure 
Maximum acceleration 

- Thermal protection system determines maximum allowable heating rates 

- Dynamic pressure and acceleration load sizes aerobrake structure 

- Heating rate and impingement angle sizes aerobrake diameter 


III. AEROBRAKE TRAJECTORY OPTIMIZATION 


A. Formulation of the Optimization Problem 

This discussion of trajectory optimization will begin with a description of the equations of 
motion which govern a trajectory and the formulation of the optimizing technique which can be 
applied to aerobrake trajectories. This will be followed by an investigation of candidate measures of 
aerobrake performance. Next, a detailed description of the process of finding the optimal angle of 
attack control function will be presented. Finally, several vehicle control strategies which could be 
used in aerobrake applications will be discussed. This section will lead to the variable angle of attack 
strategy used in section IV. 

The motion of an aerobraked vehicle can be described by seven differential equations which 
apply during both the space portion and atmospheric portion of a trajectory. These equations are 
written below in vector form [1]: 


where 


r = V, 

(3-1) 

V = (T+A)/m + g(r,t) , 

(3-2) 

m = T/goI S p , 

(3-3) 

zL 

! 

II 

00 

(3-4) 
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From these equations it can be seen that during the space portion of the trajectory the motion is 
described by basic Keplerian orbital mechanics since ITI = IAI = 0. The vectors T and A (thrust and 
aerodynamic forces) are used as controls to achieve a desired transfer from one orbit to another. The 
initial and final orbits are represented by the boundary conditions: 


A 

S 

ii 

3 

(3-5) 

V(f 0 ) = V 0 , 

(3-6) 

m(to) = "JO . 

(3-7) 

G[r (t f ), V (tf)] = 0 . 

(3-8) 


G is a vector-valued function, such as the set of orbital elements at atmospheric exit which describe 
the final orbit. The orbit transfer is described by defining the functions T(r) and A (t) during the time 
interval to<t< tf that transfer the vehicle from its initial state, equations (3-5) to (3-7), to the 
desired final state, equation (3-8) [1], 

The function T(r) represents the thrust level during propulsive bums. As the number and 
duration of the burns increase, the propellant consumption increases according to equation (3-3). An 
alternative measure of propellant consumption can be derived from equation (3-2). This measure is 
called the characteristic velocity, C, and is defined by integrating the thrust acceleration over the 
entire transfer [1]: 


T = T/m , 


(3-9) 


C = 



(3-10) 


The characteristic velocity is the sum of all of the propulsive velocity changes. It provides a more 
general measure of propellant consumption since it does not depend on specific impulse [1], For this 
reason, the characteristic velocity has been a popular performance index for evaluating aerobraked 
orbital transfers. 


It was stated earlier that the control functions T(r) and A (t) are chosen to achieve a desired 
orbital transfer. In most cases, there are many functions which will result in the same transfer. The 
analysis of aerobrake trajectories then must include an evaluation of which functions provide the best 
transfer as indicated by some performance criteria. In order to select the best control functions, an 
optimization process is required. The optimization techniques applied to aerobrake trajectories make 
use of the equations of motion and a performance index to formulate an optimal control problem. 
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The selection of the optimization formulation depends somewhat on the performance index. 
Previously, the characteristic velocity was mentioned as a possible performance index. Other 
potential performance indices include aerodynamic heating rate, acceleration load, dynamic pressure, 
or altitude drop. If the performance index is a parameter which must be integrated over the duration 
of the aeropass, it is convenient to write the variational formulation using the Bolza form of the 
functional J [1]: 

( tf 

J = S[r(t f ), \(t f ), m(tf )] + F[r(r), V(r), m(r)] dt . 

jfo (3-11) 

The function 5 represents the final boundary conditions that are possible from the allowable trajec- 
tories. The function F is the trajectory or vehicle parameter which is to be integrated over the 
aeropass. F is a function of position, velocity, and mass which in turn are functions of the general 
control function u (t). The allowable trajectories are determined from two sets of auxiliary equations 

[ 20 ]: 

xi = 0 and ki = 0. (3-12) 

The equations xi represent the governing equations for the system. In this case, these are the 
equations of motion, equations (3-1) to (3-3). The equations k t represent boundary constraints such 
as the maximum acceleration or heating rate that is permissible for the trajectory. In some 
formulations [1], this is written as an inequality constraint; 

*i[r(f), v(r)] < 0 . (3-13) 

This makes it easier to undertand the role of this function. As long as the value of A, is below the 
maximum value, Jk; plays no part in determining the optimal trajectory [1], but at the maximum value, 
ki acts as a constraint on the trajectory. 

The goal of the optimization is to determine the control functions T(r) and A(r) and the cor- 
responding trajectory which will minimize J. This is accomplished by defining a Hamiltonian in the 
following form [20]: 

H = F+'Zp i x i . (3-14) 

The pi are adjoint functions for each degree of freedom. They can be components of an adjoint vector 
or adjoint scalars. They serve the same purpose as Lagrange multipliers in an Eigenvalue problem 
except that they are functions of time rather than constants [20]. Using equations (3-1) to (3-3) and 
equation (3-11) the Hamiltonian can be written as [1]: 


H = F{r,X,m) + p r • V + p v • [(T+A)/m+g] - Pm Tlgo hp - C - 0 • (3-15) 

I I I 

(3-1) (3-2) (3-3) 
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Since time does not appear explicitly, the Hamiltonian is equal to zero [1]. In the above formulation, 
p r and p v are adjoint vectors (each with three components), and p m is a scalar. The functional J is 
minimized by minimizing the value of H at each point on the trajectory. In this case, the minimum 
value is equal to zero. The Hamiltonian is minimized by satisfying the Euler-Lagrange equations 
which relate H to the adjoint functions [20]. 


jdH 

lax, 


_d_jdH 

dx i lax, 


= 0 i- 1 ,n . 


(3-16) 


The xi represent the generalized coordinates (dependent variables) of the problem (r x , r y , r z , v x , Vy, 
v z , m). Since the entry and exit conditions for the aeropass may be variable, the transversality 
condition shown below must also be satisfied [20]: 



. dH \ . y< dH , 


+ dt = o . 


(3-17) 


Equations (3-16) and (3-17) and equations of motion y/ are solved using the control functions as 
independent variables for all of the unknowns, xj(t) and p,(r). 

If the performance index represents a parameter for which a peak value is to be minimized 
(regardless of where in the trajectory the peak occurs), the functional is written in a "nonclassical" 
form as a "minimax" or "Chebyshev" optimal control problem [21]: 


J = S[r(tf), m(tf )] + max F[r(r), V(r), m(t )] . 


(3-18) 


To simplify the notation, this can be rewritten as: 


7(u) = S(\(t/)) + max F(x,u) , (3-19) 


where u is the control function (u replaces T and A), and x is the set of state vectors. Using this 
notation, the equations of motion can be abbreviated by: 

X = f(x,u) . (3-20) 


The functional is minimized by satisfying a number of necessary conditions during the trajectory. The 
Hamiltonian is defined as [22]: 


//(x,p,u) = p r f(x,u) . 


(3-21) 
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The p represents adjoint functions similar to those used in the Bolza form. 

The necessary conditions for an optimal trajectory can be described by a set of equations that 
must be satisfied at different times during the aeropass depending on whether the performance index 
is at its peak value. One condition that must be satisfied at all points on the trajectory regardless of 
the value of the performance index is the condition for the Hamiltonian [22]: 

tf(x,p,u) = c , (3-22) 


(if time is not specified explicitly, H = 0). 

If the performance index, F, is below its peak value, the following condition must be satisfied 

[ 22 ]: 


P (f) = - 


BH_ 

dx 


(3-23) 


At points in the trajectory where F hits peak values (assuming it is possible to have multiple peaks 
at t\, t2,...tk), another equation must be satisfied [22]: 


dF 

p(r + ) = p(r ) + p, — (xfo)). 


i = 1,2,...,* , 


(3-24) 


where 

k 

X P= 1 . p>0. 

i=l 


This condition relates the value of p just prior to the peak to the value just after the peak, and it is 
used to account for the discontinuity in p(r) at the peaks. If the value of F remains at the peak value 
for some finite time interval (//, t\), the trajectory has a "flat maximum" [22], During this time 
interval two other equations must be satisfied [22]: 


• (t , _ dH . R 3^(x,u) 

P (0 - - + o, — ^ » 

dx dx (3-25) 


where 8 ,(r) is a scalar function satisfying: 

dH s . . d F 
T-=5iWT-. 

du du (3-26) 
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Finally, there is a necessary condition that must be satisfied at atmospheric exit [22]: 


P (tf) = - a 



(3-27) 


where ot is a set of multipliers. A summary of which equations must be satisfied at various times to 
obtain an optimal trajectory using the Chebyshev formulation is shown in table 3-1. 


Table 3-1. Necessary conditions for a Chebyshev optimal control problem. 


PERFORMANCE INDEX VALUE EQUATIONS TO SATISFY 

F< MAXIMUM VALUE (3-22), (3-23) 


F REACHES MAXIMUM (3-22), (3-24) 


F IS AT A FLAT MAXIMUM (3-22), (3-25), (3-26) 

AT ATMOSPHERIC EXIT (3-22), (3-27) 


B. Performance Indices 

The optimum trajectories that are found using the formulations that have been described pro- 
vide the best performance as defined by the selected performance index [1], Recent work has 
examined a number of performance indices and the implications they have on the resulting optimal 
trajectories [3]. It is possible that two performance indices will yield completely different "optimal" 
trajectories. A list of possible aerobrake trajectory performance indices is shown in table 3-2. The 
choice of a performance index often depends on specific mission or vehicle requirements An aero- 
brake that is designed to perform a plane change usually flies deeper into the atmosphere than an 
aerobrake designed only to decrease velocity. Since the trajectories for these two purposes are 
inherently different, it would probably be necessary to use different performance indices to optimize 
each type of mission. The choice of a performance index must be made carefully. It is possible that if 
a trajectory is optimized to satisfy one criteria, the vehicle may be penalized in another aspect. The 
choice of a performance index often involves a compromise of several vehicle and trajectory 
constraints. J 


The most obvious performance index for an aerobrake trajectory is the characteristic velocity, 
C. It was defined earlier by equation (3-10). An alternative definition replaces the integral in equa- 
tion (3-10) with a simple summation of the delta velocities [3]: 

C ~ A Vjotai = A V\ +AF 2 +A V , '}+... . (3-28 

In the optimization formulation, equations (3-10) or (3-28) are substituted for F into the functional. 
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Table 3-2. Aerobrake performance indices. 

1 . Characteristic velocity (or propellant mass) 

2. Payload mass fraction 

3. Integrated convective heating rate 

4. Peak convective heating rate 

5. Peak acceleration 

6. Peak dynamic pressure 

7. Time integral of flight path angle squared 

8. Peak altitude drop 


Minimizing the characteristic velocity is in many cases equivalent to minimizing the propel- 
lant usage [1] or vehicle mass. This implies the minimum mission cost since cost is often propor- 
tional to vehicle mass [1]. Since the characteristic velocity is primarily a measure of propellant 
requirements, it overlooks two other factors that comprise the vehicle mass: the payload mass and 
the vehicle inert mass. The payload is the cargo or instrument that is carried by the veicle. The inert 
mass consists of the vehicle structure, subsystems, and aerobrake. Using the above three factors, 
the total vehicle mass can be written as [1]: 


^vehicle = ^propellant + m pay load + m inert • 


(3-29) 


If an aerobrake mission is to be advantageous, the decrease in propellant mass must be greater than 
the increase in inert mass due to the addition of the aerobrake [1]. 

The balance between propellant savings and aerobrake mass is more obvious if another per- 
formance index is used. The payload mass fraction sometimes provides a better measure of 
aerobrake performance than the characteristic velocity [1]. It is defined as: 


M — Wpayioad/ ^vehicle ■ 


(3-30) 


This ratio automatically allows for the trade-off between propellant mass and aerobrake mass. The 
oDtimum trajectory, defined by the maximum value of M, would provide the best combination ol 
energy requirements (delta velocity) and aerobrake requirements (TPS). The payload ratio implies 
that there is an inherent relationship between the vehicle design (specifically aerobrake TPS mass) 
and the trajectory design. 


The coupling between TPS mass and trajectory design has been considered in many studies 
[1]. Performance indices which characterize this relationship more precisely than payload ratio have 
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been investigated. In most cases, the relationship is not direct but in terms of some measurable 
trajectory parameter. Some of these alternative performance indices are listed in table 3-2. The peak 
heating rate (0max) and the total heat load ( Q ) are the most common indices used as an indication of 
TPS system requirements. 

The peak heating rate is computed at the vehicle stagnation point [1] and is defined by the 
relation [1]: 


Q = c p 0 5 V 3 08 W/cm 2 , 


(3-31) 


where 


c= 1.83 x 10- 8 r n -0.5 (l-g w ) 
r n = nose radius (m) 
g w = wall enthalpy/total enthalpy 
p = density (kg/m 3 ) 

V = velocity (m/s) . 


In many cases c is dropped from the functional, so the performance index has the form [3]: 


F = p0.5y 3.08 


(3-32) 


The integrated heat load is simply the heating rate integrated over the duration of the aeropass [1]: 

0 .rv 

-wo (3-33) 

Two other performance indices which serve as indicators of vehicle structural requirements are the 
vehicle acceleration (V) and the dynamic pressure. The dynamic pressure gives a measure of 
aerodynamic loads. It is defined as: 


9 = 



(3-34) 


The four performance indices just described give a more precise indication of aerobrake or vehicle 
requirements than the payload ratio, but they do not necessarily lead to the best "overall" com- 
promise between vehicle and trajectory design factors. The last two performance indices listed in 
table 3-2 were proposed as indices which could be used to obtain optimal trajectories in which the 
conflicting aerobrake mission and vehicle requirements could be effectively balanced. 
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The time integral of the flight path angle squared was proposed as a performance index by 
Miele [1]. It leads to a "nearly grazing trajectory" which will be discussed later. This performance 
index is defined by the integral [3]: 


F = I xy 2 dt , 

ho (3-35) 

where x is the normalized time in the atmosphere (= 0 at entry, = 1 at exit). The nearly grazing 
trajectory seeks to barely skim the top of the atmosphere. Miele has shown that this type of trajec- 
tory has several advantages [3]. Vinh and Lu [22] proposed the minimum peak altitude drop as a 
performance index which achieves the same goal as the integrated square of the flight path angle. 

The peak altitude drop refers to the depth of penetration into the atmosphere. This performance index 
yields trajectories that are similar to the "nearly grazing" solutions. 

Miele and his colleagues [1] have conducted an extensive investigation of a number of per- 
formance indices in order to determine which leads to the most favorable flight path along with 
having favorable propellant requirements for the overall orbit transfer [1]. The indices that were 
evaluated are listed below [1]. 

1. Characteristic velocity (or propellant mass) 

2. Integrated heating rate 

3. Peak heating rate 

4. Time integral of the flight path angle squared. 

The trajectory that results in the minimum characteristic velocity is called the grazing trajectory [1]. 
It is an idealized trajectory which is not really possible, in which the vehicle travels along the edge of 
the atmosphere until just enough energy is dissipated to reach the final target orbit (fig. 3-1). The 
idealized "grazing" trajectory gives a good approximation for propellant requirements if the initial 
orbit is high or the entry angle is small and the exit flight path angle is near zero [1]. This type of 
trajectory is useful in many aerobrake performance analyses and will be discussed further in section 
IV. 


The second performance index, the integrated heating rate, is minimized by penetrating 
deeper into the atmosphere and dissipating the energy as quickly as possible [1]. This is opposite to 
the strategy used to minimize the characteristic velocity (grazing trajectory). However, even in this 
case, where the deepest penetration into the atmosphere is desired, the propellant requirement is 
only 1 percent greater [1]. The disadvantage of this performance index is that it leads to optimal 
trajectories with very high peak heating rates. This implies greater TPS requirements, thus minimiz- 
ing the "integrated" heating rate, and the "peak" heating rate yield conflicting requirements. The 
proper compromise depends on the type of TPS [1]. 

The third index, peak heating rate, is a function of density (altitude) and speed. The speed 
goes through the same range of values for any allowable trajectory, therefore the peak heating rate is 
minimized by minimizing the density during the aeropass [1]. This can be accomplished by flying as 
close to the atmospheric boundary as possible. As a result, the idealized or "grazing" trajectory 
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Figure 3-1. Idealized aeroassisted transfer. 

minimizes both the propellant consumption and the peak heating rate [1], Unfortunately, it is 
impossible to fly a grazing trajectory if realistic control bounds are applied to the optimization. The 
grazing trajectory is useful because it provides an ideal or goal for which actual trajectories can be 
designed to approximate. To achieve the best possible approximation, another performance index 
was needed. 

The desire to approximate the idealized optimal trajectory led Miele to formulate the fourth 
performance index and to name the resulting optimal trajectory the "nearly grazing trajectory" [1], 

By minimizing the flight path angle at each point in the trajectory, the penetration into the 
atmosphere is minimized. The actual altitude profile for the nearly grazing trajectory is determined by 
the trajectory boundary conditions and the vehicle control capabilities. A comparison of several tra- 
jectory parameters for transfers from GEO to LEO that result from the different performance indices 
is shown in figure 3-2 [3J. 

An interesting result of Miele's work is that performance indices 1, 3, and 4 yield almost 
identical trajectories and trajectory control functions. Miele expresses this result with the statement, 
"In the coplanar case, what is good for peak heat rate is good for the energy." From the data shown 
in figure 3-2, it appears that the nearly grazing solution provides an excellent compromise between 
energy requirements and aerodynamic heating requirements for GEO to LEO transfers [3]. The 
nearly grazing solution requires about 3.1 percent more energy than the minimum energy trajectory 
but reduces peak heating by 27.1 percent, dynamic pressure by 42 percent, and peak acceleration by 
31.6 percent [3]. Miele also has shown that the nearly grazing solution is the best compromise for 
other types of transfers as well [3]. 
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CHARACTERISTIC INTEGRATED FLIGHT PEAK HEATING 

VELOCITY PATH ANGLE RATE 

PERFORMANCE INDEX 

Figure 3-2* Relationship between performance index and trajectory parameters* 

Vinh and Lu have arrived at results similar to Miele using the minimum peak altitude drop as 
a performance index [22]. The minimum peak altitude drop is defined as the minimum altitude during 
the aeropass (periapsis). In their analysis, Vinh and Lu use the Chebyshev minimax problem to 
prove that the heating rate and dynamic pressure pass through their respective maximum values 
along the descending portion of the aeropass, and that these parameters can be reduced by 
minimizing the peak altitude drop. They also derive the necessary control function u(t) needed to 
achieve the optimal (nearly grazing) trajectory. The near optimal trajectories that will be analyzed in 
the following section will be based on the results of Vinh and Lu, therefore an overview of then- 
analysis will be presented. 


C. Lift Control Function for Optimum Trajectories 

The first step of the analysis was to describe the motion of the aerobraked vehicle in terms of 
dimensionless equations of motion. The use of modified Chapman's variables allows the problem to 
be formulated such that the only vehicle parameters are the maximum L/D and maximum lift 
coefficient. The modified Chapman's variables are defined as [23]: 

__ PSQ* JT 

Z ~ 2m V B (3-36) 
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( 3 - 37 ) 


v = 


vl 

& ’ 


5 



y cos y dt . 


( 3 - 38 ) 


The atmospheric density is assumed to be inversely proportional to the scale height B, 


dp = -Bp dr . 


( 3 - 39 ) 


The variable z is the altitude drop. It increases as altitude decreases and can be regarded as a 
replacement for altitude. The variable v is the dimensionless kinetic energy which is a measure of 
speed. Finally, the variable s is the dimensionless arc length of the atmospheric trajectory. It is the 
independent variable [22], 

The vehicle parameters that will be used in the nondimensionalized equations of motion are 
derived from an assumed parabolic drag polar for the vehicle [22]: 


c d = c Do + kc l 2, 


at the maximum LTD, E *: 


( 3 - 40 ) 


therefore: 


= and C D * = 2 C Do , 


= — k , 

2 UcCoq 


( 3 - 41 ) 


( 3 - 42 ) 


the normalized lift coefficient, which is the control and is defined as: 


*•=0? • c * sc '™- 


( 3 - 43 ) 


Using the above definitions, the dimensionless equations of motion for the nonthrusting 
vehicle considered as a point mass were written as [22]: 
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(3-44) 


= -k 2 z tan y , 
ds 


dv_ _ -*zv(l+^ 2 ) _ , 2 _ 


ds E* cos y 


(2-v) tan y , 


(3-45) 


_ kz X cos q + ( i _ 1\ 
ds cos y ' vP 


(3-46) 


dQ = cos ¥ 

ds cos <(> ’ (3-47) 


d<j> 

— = sin 
ds 


(3-48) 


_ /c; 3, sin C L _ cos y tan <[> , 

ds cos 2 Y (3-49) 


where: 

Y = flight path angle 
0 = longitude 
<p = heading angle 
X = normalized lift coefficient 
c = bank angle. 

The quantity k : 2 = Br is treated as a constant with a value of 900 for Earths atmosphere [22]. 

Vinh and Lu chose three trajectory parameters as possible performance indices. The first was 
the stagnation point heating rate, defined as: 


Q = p0.5 V 3.08 


(3-50) 
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The second performance index was the dynamic pressure, defined as: 


q = 



(3-51) 


The third performance index was the peak altitude drop, z, which was defined earlier by equation 
(3-36). By substituting the Chapman variables into equations (3-50) and (3-51) and rescaling the 
formulas, all three performance indices can be put into the form [22]: 


F = K zv n , 


(3-52) 


where k is a constant. If n = 3.08, F represents peak heating rate; if n = 1, F represents dynamic 
pressure; and if n = 0, F represents peak altitude drop. 

Using the derivative of F with respect to time, it is possible to determine that the maximum 
value of F for any value of n,n> 0, occurs during the descending part of the trajectory. 


dL - dF_ ds, _ fr v „-i Ldz + nz dv\ds_ 

dt ds dt ' ds ds' dt' (3-53) 


from the definition of s, (equation (3-38)) 


^> 0 , 

dt 


also, at atmospheric entry: 


z = 0 , 


<&> 0 , 

ds 


(z increases as altitude decreases). Due to the conditions of equations (3-54) to (3-56), it can be 
concluded that at atmospheric entry [22]: 

^> 0 . 

dt (3-57) 


At the lowest point of the trajectory, since the flight path angle is zero: 


(3-54) 

(3-55) 

(3-56) 
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<& = 0. 
ds 


(3-58) 


Since the speed is always decreasing: 


^< 0 , 

ds 

therefore at periapsis, equation (3-53) will result in: 

^< 0 . 

dt 


(3-59) 


(3-60) 


Since the time derivative of F is positive at atmospheric entry and negative at periapsis, it must 
pass through zero somewhere on the descending path. This point corresponds to a maximum value of 
F. From equation (3-53), it can also be concluded that the location of the maximum value of F 
depends on the value of the exponent n [22]. Vinh and Lu have proven that the peak heating rate 
always occurs first, followed by the peak dynamic pressure, then the peak altitude drop [22], From 
the trajectory properties just described, it can be inferred that by minimizing the peak altitude drop, 
the peak heating rate and peak dynamic pressure are also reduced [22]. 

The process of finding the control function which yields the minimum peak altitude drop begins 
with the formulation of the optimization problem. If only coplanar aeroassisted transfers are con- 
sidered, the dimensionless equations of motion, equations (3-44) to (3-49), reduce to [22], 


dz = - it 2 z tan y , 

ds 


(3-61) 


dv _ ~ kzv _ (2-v) tan 

ds E* cos X 


Y. 


(3-62) 



(3-63) 


The aerodynamic control is the normalized lift coefficient defined by equation (3-43). The perform- 
ance index in its most general form is given by equation (3-52). The Hamiltonian is defined using 
equations (3-61) to (3-63) [22]: 


H = —pi k z tan y - 


+ <2-v) “ Y ) +ft (^v + M)= 0 • 


(3-64) 
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Subject to the necessary conditions that are outlined in equations (3-22) to (3-27), the 
optimal lift control can be determined. For an isolated maximum, the lift control is defined by the 
following relation [22]: 


X = 


E *Py 
2 vp v 


, if |X| < X 


^ma x sign (py ) , 


At each point in the trajectory, the adjoint variables 


max 

otherwise 


must satisfy the following equations [22]: 


(3-65) 


dH ,2 

Pz = --z- = Pz^ tan y + 

az 


p v kv{ 1+X 2 ) PykX 
E* cos y ~ cos y ’ 


(3-66) 


dH_ = p 1 kz{ 1+X 2 ) 
dv E * cos y 


p v tan y - 


Py 


(3-67) 


dtf I 2 Pv kzv (1+X 2 ) sin y] i . , , , i 

Py = ~^Z = \Pz K z + + (p v (2-v) -py kz X sin y) — 1 

' E* 'cos 2 y ' cos 


cos 2 y 


(3-68) 


The optimal trajectory is defined by integrating the system of equations (3-58) to (3-60) and (3-66) 
to (3-68) using equation (3-65) [22], 

If the performance index has a flat maximum, then during the time at the maximum equation 
(3-65) is replaced by [22]: 


^2 _ -E* sin y (m(2-v) + k : 2 v) ^ 
m kzv 


(3-69) 


and the adjoint variables must satisfy the following equations instead of equations (3-66) to (3-68) 

[ 22 ]: 


Pz = 


dH ~dF 

-3T + 8 37' 


(3-70) 
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(3-71) 


Pv = 


Py ~ 


5 = 


p - dE_ _ 

ds 


dH 

i 


dv 

b dv ’ 

dH 

5 — 



dy 

dy ’ 

Py V X 

~E*Py 

2 n k z 

V"TX ’ 

nkzv 

(1+A 2 ) 

\ E* 

cos y 


(3-72) 


(3-73) 


(3-74) 


It was mentioned earlier that by minimizing the peak altitude drop it is possible to reduce 
both the peak heating rate and the peak dynamic pressure. This assertion enables the formulation to 
be simplified by allowing the value of n to be set to zero to obtain the best overall optimal trajectory. 
The performance index is simplified to [22]: 


F = kz. (3-75) 


The use of peak altitude drop as the performance index also significantly simplifies the Hamiltonian 
and adjoint functions. Vinh and Lu have applied numerical methods to solving the Chebyshev prob- 
lem that has just been described, and they have defined the optimal lift control functions for several 
sets of atmospheric entry and exit conditions using the peak altitude drop as the performance index. 
The lift control function that they have defined for coplanar aeroassisted transfers from high altitude 
orbits to low altitude orbits will be the basis for the angle of attack profiles used for the trajectories 
in section IV. 

The lift control profile and the resulting altitude profile for the optimal trajectory defined by 
Vinh and Lu are shown in figures 3-3 and 3-4. This specific trajectory was calculated using the 
following values for the maximum normalized lift coefficient and maximum L/D [22]: 

X = 2.0 and E* = 0.5 . 

At atmospheric entry, the lift coefficient is at its maximum positive value (lift up). It remains at this 
value until shortly before the peak altitude drop (periapsis). This allows the vehicle to fly at the 
highest altitude possible during the descending portion of the trajectory. Since the peak heating rate 
and peak dynamic pressure occur during the descending leg of the trajectory, it is beneficial to fly at 
the highest possible altitude where the density is less. Shortly before periapsis, the lift control 
begins a continuous transition toward the maximum negative value of lift coefficient (lift down). The 
lift coefficient has a value of zero at periapsis [22], During the ascending leg of the trajectory, the lilt 
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Figure 3-3. Variation of lift control for a minimax aerobrake trajectory [22]. 
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Figure 3-4. Altitude drop for minimax trajectory [22], 
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control remains negative but decreases in magnitude. The negative lift allows the vehicle to remain 
in the atmosphere as long as possible and also have a small flight path angle at atmospheric exit. 
The decrease in magnitude of the lift coefficient results in an increase in the drag coefficient so the 
maximum velocity change can be achieved. 


D. Control Strategies 

Based on the results of Miele, Vinh, and Lu, it appears that the nearly grazing solution is the 
best "overall" optimum aerobrake trajectory. The aerodynamic control function can be obtained using 
the integral of the flight path angle squared [1] or the minimum peak altitude drop [22] as the per- 
formance index in the optimization formulation. The aerodynamic control functions are usually speci- 
fied by the magnitude and direction of the lift force and the magnitude of the drag force. In order to 
describe these functions in terms of vehicle characteristics and attitude conditions, it is necessary to 
choose a vehicle control strategy. There are three possible control strategies that can be used for 
aerobrakes. They are roll control, pitch control, and drag control (sometimes called drag modulation). 
In theory, they could be used separately or in combination). 

Roll control uses a rotation of the vehicle about its roll-axis (bank angle) to control the direc- 
tion of the lift vector. By rotating the lift vector out of the orbital plane, the vertical component of lift 
is used to control the rate of altitude change, and the lateral component of lift is used to control the 
rate of heading change (plane change). The bank angle profile for an aeropass for a GEO to LEO 
transfer is shown in figure 3-5 [3]. This profile shows the bank angle is varied and below 90°. This 
results in a 22.8° plane change in the atmosphere. If no plane change is desired, the bank angle is 
alternated from side to side such that the time spent banked to each side is balanced so there is no 
net heading change. One method for doing this is to perform a series of S-tums during the aeropass. 
Another method is to continuously roll the vehicle at varying rates during the aeropass [24], The two 
strategies for roll control are illustrated in figure 3-6. In both strategies, the rotation of the vehicle 
would be achieved using the reaction control system (RCS) thrusters on the vehicle. 

Lift 



Normalized Time of Flight x 

Figure 3-5. Bank angle profile for a GEO to LEO transfer [3]. 
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Figure 3-6. Roll control strategies. 
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Roll control requires a very responsive control system. The RCS system will need roll accele- 
ration rates of about 107s 2 [14]. The thrust required to achieve this level of angular acceleration is 
shown in figure 3-7. The thrust level depends on the vehicle inertia in the roll axis and the moment 
arm of the thrusters. The propellant required by the RCS system is shown in figure 3-8. The data 
shown is for a return from lunar orbit and assumes 12 roll maneuvers at an L/D of 0.15 during the 
aeropass [14]. 

Pitch control uses a rotation of the vehicle about its pitch axis (angle of attack) to control the 
magnitude of the lift and drag forces. There are no lateral forces associated with pitch control, there- 
fore, pitch control cannot be used for atmospheric plane changes unless it is used in combination with 
roll control. The angle of attack for an aerobraked vehicle could be controlled by mechanically shifting 
the vehicle center of mass such that the vehicle would "trim" at the desired angle of attack. Two 
strategies for changing the vehicle center of mass are shown in figure 3-9. One strategy would use 
tracks or sliding mechanisms to shift the position of certain vehicle elements. The other strategy 
proposed by Tom Goodrick of Marshall Space Flight Center (MSFC) is to rotate the aerobrake pay- 
load with respect to the aerobrake using a pivot mechanism. This strategy has the advantage that 
the payload is kept further from the flow behind the aerobrake. 

Pitch control may offer several advantages over roll control. It allows a continuous, smooth 
variation of the L/D from 0 to as much as 1.5, if flow impingement can be managed [13]. Another 
advantage is that all of the motion remains in the orbit plane, so many of the calculations are simpli- 
fied [13]. Since the vehicle does not require S-tums or a continuous roll, the vehicle motion is 
reduced. This may simplify the guidance scheme and reduce the requirements on the RCS system. 

The third aerobrake control strategy is drag control. It uses a change in the surface area or 
shape of the aerobrake to control the magnitude of the lift and drag forces. As was the case with 
pitch control, there are no lateral forces generated so drag control cannot be used for atmospheric 
plane changes unless roll control is also used. The primary difference in control authority between 
pitch control and drag control can be seen in the equations for lift and drag: 

L=J-p V 2 C l A D = lp V 2 C d A. 

2 V 2 K (3-76) 


Pitch control uses angle of attack to vary the lift and drag coefficients. Drag control changes the 
vehicle reference area without necessarily changing the aerodynamic coefficients. The effectiveness 
of one strategy compared to the other depends on the relationship between the aerodynamic coeffi- 
cients and angle of attack and the feasibility of changing the aerobrake surface area. One concept 
that has been proposed for controlling the aerobrake area is to use an inflatable elastic aerobrake 
(called a ballute). An illustration of a ballute aerobrake is shown in figure 3-10 [25]. 

Each of the three control strategies has advantages and disadvantages. The ideal control 
strategy would probably make use of all three types of control, but when realistic vehicle consider- 
ations are imposed this may not be possible. There has been considerable study of aerobrake tra- 
jectories which use only roll control. One possible reason is that roll control has been used almost 
exclusively for atmospheric reentry vehicles dating back to Project Mercury. Many of the computer 
programs used for reentry analysis can be modified to study aerobraking to orbit. A possible reason 
that pitch control has not received a high level of study in recent years could be the fact that the 
planned AFE will use only roll control. Pitch control also has the design problems posed by the 
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Figure 3-10. Ballute aerobrake [25]. 

control mechanisms which have hindered its popularity. Drag control was studied during the early 
1980's in NASA’s Space Transfer Vehicle studies, but the technological problems of the ballute have 
diminished its popularity. 

Despite the fact that pitch control has not been as popular as roll control in recent aerobrake 
studies, there are several reasons that it may be more effective. These reasons can be illustrated by 
examining the differences between aerobraking to the surface (reentry) and aerobraking to orbit. 
These differences are summarized in figure 3-11. The most obvious difference is in the trajectory 
paths. In aerobraking to the surface, the trajectory becomes steeper during the reentry, while in 
aerobraking to orbit the trajectory remains almost parallel to the surface. Probably the most signifi- 
cant difference is in the velocities. In aerobraking to the surface, the velocity quickly decreases to 
suborbital speed, but in aerobraking to orbit, the velocity is always above orbital velocity. The target 
condition for a reentry trajectory is usually a specified landing site. Since the trajectory is always 
becoming more vertical, the reentry trajectory becomes less sensitive to velocity errors as long as 
the vehicle maintains the proper path. The target condition for an aerobraking to orbit is usually a 
specified apoapsis altitude at atmospheric exit. This altitude is very sensitive to both exit velocity 
and exit flight path angle. A small error in velocity due to guidance errors or atmospheric density 
variations results in a greatly amplified error in final apoapsis altitude. 
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Figure 3-11. Differences between aerobraking to the surface and aerobraking to orbit. 

The differences in sensitivity to velocity errors between reentry trajectories and aerobraking 
to orbit trajectories indicates the possibility that the two different types of aerobraking may require 
different guidance and control strategies. Trajectories for aerobraking to orbit appear to be much 
more sensitive to velocity errors. The control system needs fast response times. If roll control is 
used, the requirement to balance the lateral components of lift will add considerable complexity to 
the guidance strategy and increase the response time to compensate for trajectory errors. The 
factors that will have the greatest influence on the response time of a control strategy will be the 
aerodynamic properties of the aerobrake and the capabilities of the control mechanism (mechanical or 

RCS). 


In the next section, the requirements of the aerobrake control system will be investigated. 
Several aerobrake trajectories will be modeled which will use variable angle of attack control strate- 
gies which approximate the optimum lift profile that has been described in this section. 


50 



IV. AEROBRAKE TRAJECTORY ANALYSIS 


A. Trajectory Simulation (IMP Computer Program) 

To investigate the relationships between aerobrake design factors and to develop approxima- 
tions of the flight characteristics of aerobraked vehicles, it is necessary to mathematically model 
their trajectories through the atmosphere. The particular methods used to perform this modeling 
depend on the type of data that is desired and the level of accuracy needed. If the only data that is 
required is a reasonable approximation of propellant consumption, detailed trajectory data is not 
necessary, but, if vehicle loads or control requirements are being investigated, then an accurate 
trajectory must be calculated. This requires careful modeling of the forces on the vehicle and the 
atmosphere in which the vehicle flies. 

The trajectory data that are included in this report were developed using numerical simula- 
tions of the aeropasses. The desired output consisted of vehicle loads and orbit conditions at 
atmospheric exit. Vehicle characteristics, entry conditions, and atmospheric density were varied in 
order to investigate the sensitivities of the trajectory and vehicle loads. The trajectory data that 
were developed will be described in detail later. 

The computer program used to model the aerobrake trajectories is called IMP (Integrated 
Mission Program). It is a simulation language used to model orbital missions about Earth, Mars, or 
the Moon [26]. It was developed by Vince Dauro of MSFC and is used in many preliminary design 
studies. IMP uses a Fehlberg seventh-order, 13-evaluation Runge-Kutta integrator with error and 
step size control [27] to integrate the equations of motion of an orbital vehicle. The equations of 
motion include the accelerations listed below [26]. 

Earth: Earth’s gravity 

Sun's gravity 
Moon's gravity 

Atmosphere (aerodynamic forces) 

Solar pressure 

Mars: Mars' gravity 

Sun's gravity 
Atmosphere 
Solar pressure. 

For aerobraking trajectories the acceleration due to solar pressure is not applicable. 

The orbital parameters are computed as smoothed mean values using up to a sixth-order 
gravity model to simulate the Earth or Mars gravity field. Earth's gravitational potential field is 
computed with J through J4 constants. The Mars gravitational field is computed using the J and J2 
constants [26]. 
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B. Atmosphere Modeling 


1. Atmosphere Models 

In order to compute the atmospheric trajectories, it is necessary to use an atmosphere model 
in the program. There are several atmosphere models which can be used within the IMP program. 
These models provide a density profile as a function of altitude. Accurate modeling of the atmos- 
pheric density is of paramount importance in aerobrake analysis since the drag is proportional to the 
density [28]. The atmosphere models which can be selected for use by the program are listed below: 

Earth: U.S. 1962 Standard Atmosphere 

Jacchia 

Mars: COSPAR Mars cool-low (temperature-pressure) 

COSPAR Mars warm-high 
MSFC Mars Gram. 

In each of the atmosphere models, the density altitudes are calculated spherically [26], In 
this report, the U.S. 1962 Standard Atmosphere was used for Earth aerobraking missions, and the 
Mars cool-low atmosphere was used for the Mars aerobraking missions. The U.S. 1962 atmosphere 
is defined in terms of an ideal air, assumed to be devoid of moisture, water vapor, and dust [29]. The 
density profile shown in figure 4- 1 was derived from the calculated values of pressure and molecu- 
lar-scale temperature. The top of the Earth's atmosphere was assumed to be 121 -km altitude. 

Above this altitude, the atmospheric drag force can be neglected for aerobrake missions. 

The COSPAR Mars atmosphere models provide envelopes around the excursions of tem- 
perature and density measured by the Viking probes during aerobraking and descent. They represent 
the best knowledge of the range of temperature and density from near 100-km altitude down to the 
surface. The density of the upper atmosphere of Mars was calculated using an estimate of the 
ultraviolet radiation on the atmosphere since the ultraviolet insolation has a greater effect on the 
exospheric temperature and density than the surface pressure [28]. The density profile is shown in 
figure 4-2, which was constructed from data in reference 30. 


2. Density Variations 

The atmosphere models that were chosen provide nominal or baseline density profiles that 
can be used in determining the sensitivity of the trajectory to vehicle characteristics or entry 
conditions. The actual atmospheric densities vary widely due to factors such as time of day or year 
and solar activity. Figure 4-3 shows that density variations of up to ±50 percent from the nominal 
value are possible in the range of altitudes where aerobraking would take place in the Earth’s 
atmosphere [2]. There is a similar range of atmospheric density variation in the Mars atmosphere. 
Aerobraked vehicles must be designed to operate over these wide ranges of atmospheric conditions 

The amount of departure from the nominal density profile is not constant with altitude. Data 
from space shuttle flights have indicated "potholes in the sky" as shown in figure 4-4. These dia- 
grams show the amount of departure from the U.S. 1962 standard atmosphere density profile. The 
Viking entry measurements showed the same type of oscillatory density variations in the Mars 
atmosphere [28]. These variations are also shown in fieure 4-4. These large density excursions, 


52 





100 




DEPARTURE FROM U.S. STANDARD DENSITY, PERCENT 

Figure 4-3. Range of Earth atmosphere density variations [2] 
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Figure 4-4. Atmospheric density variations during reentry. 

which occur over relatively narrow altitude ranges, pose a significant challenge in the design of aero- 
brake control systems [2], Since it will be difficult to predict the atmospheric conditions ahead of the 
vehicle, sophisticated adaptive control laws may be required [2]. 

Recently, aerobrake guidance and control system design has received intensive study. The 
central issue in guidance and control system design is the amount of aerodynamic control authority 
required to contend with unpredictable variations in atmospheric density. This report will not attempt 
to define a guidance law that could be used in an actual aerobrake mission. Instead, only the bounds 
of the requirements for the control systems will be investigated. This will be accomplished by 
investigating the capabilities of low L/D aerobrakes for several aerobrake applicationns. As an initial 
step in this analysis, it is necessary to determine the sensitivity of aerobrake trajectories to entry 
conditions, vehicle characteristics, and atmospheric conditions. The following section will show tra- 
jectory data that were generated to investigate the relationships between several aerobrake 
parameters. 
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C. Aerobraking Trajectory Analysis 

Numerical simulations of aerobrake trajectories were used to analyze four applications of 
aerobraking. These applications were return from geosynchronous orbit, return from lunar orbit. Earth 
aerocapture from a hyperbolic orbit, and Mars aerocapture from a hyperbolic orbit. The analysis 
consisted of identifying the sensitivity of the apoapsis altitude at atmospheric exit to the incoming 
target vacuum periapsis altitude and the vehicle angle of attack. The vacuum periapsis represents 
the periapsis altitude of the orbit if there was not an atmosphere present. The peak vehicle loads 
during the aeropass were also investigated. These loads included the peak acceleration load, the 
peak stagnation point temperature, the peak dynamic pressure, and the peak stagnation point con- 
vective heating rate. The sensitivity of these loads to entry conditions, ballistic coefficient, and angle 
of attack was investigated using a wide range of constant and variable angle of attack trajectories. 
The analysis was completed with an investigation of the effect of atmospheric density variations on 
the range of allowable target periapsis altitudes and the definition of a target vacuum periapsis 
tolerance to allow for these variations. 

The data for return from geosynchronous orbit are representative of atmospheric entry from an 
elliptical transfer orbit with an apoapsis altitude of 35,876 km and a target vacuum periapsis inside 
the Earth’s atmosphere. The data for return from lunar orbit is representative of atmospheric entry 
from a transfer orbit from the Moon, with an apoapsis altitude of 385,747 km. The data for Earth 
aerocapture covers the range of entry conditions for an atmospheric entry from a hyperbolic orbit 
representative of a return from Mars. Similarly, the Mars aerocapture data covers the range of entry 
conditions representative of a transfer from Earth. 

Table 4-1 lists the orbital and vehicle characteristics that were used in the trajectory simu- 
lations. The range of ballistic coefficients and mission energies (represented by C 3 s) were derived 
from studies performed at MSFC from 1983 to the present. The aerodynamic properties of the vehicle 
were based on the spherical, Apollo style aerobrake. The angle of attack range yields low L/D s 
between 0.0 and 0.25, with the lift directed upward or downward. The heating rate calculations are 
based on the stagnation point convective heating rate defined by equation (3-31). A nose radius of 
15 m was used. 


1 Sensitivity of Trajectory to Entry Con d itions and Vehicle Characteristics 

Figures 4-5 through 4-11 show the sensitivity of the orbit apoapsis altitude at atmospheric 
exit to target vacuum periapsis altitude for various ballistic coefficients and angles of attack for each 
of the aerobrake applications. These data were generated assuming a constant angle of attack for the 
entire aeropass. In reality, this would not be possible for reasons that will soon become obvious. The 
data for angles of attack of -15° (lift upward) clearly show that as the target vacuum periapsis in- 
creases, the exit apoapsis increases. Also, the exit apoapsis becomes much more sensitive to target 
periapsis as the periapsis altitude increases. It approaches a point where the exit apoapsis is so 
large that it is difficult to see a relation between target periapsis and exit apoapsis. This sensitivity 
increases as angle of attack becomes more positive. In most cases, the curves for angles of attack of 
10 ° or 15° (lift down) are almost vertical lines. 

The extreme sensitivity of exit apoapsis to target periapsis altitude indicates that it would 
not be realistic to use a constant angle of attack (constant L/D) trajectory. This is especially true in 
the case of aerobrake applications where the desired exit apoapsis altitude is high. In many aero- 
capture studies, it is assumed that the vehicle would aerobrake into a highly elliptical orbit with 
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Figure 4-6. Lunar return aerobraking — exit apoapsis as a function of target vacuum periapsis 
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Figure 4-7. Earth aerocapture — exit apoapsis as a function of target vacuum periapsis 

C 3 = 10 km 2 /s 2 . 
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Figure 4-8. Earth aerocapture exit apoapsis as a function of target vacuum periapsis at 

C3 = 50 km 2 /s 2 . 
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Figure 4-10. Mars aerocapture — exit apoapsis as a function of target vacuum periapsis at 

C 3 = 25 km 2 /s 2 . 
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Figure 4-11. Mars aerocapture — exit apoapsis as a function of target vacuum periapsis 
° C 3 = 40 km 2 /s 2 . 




a period of 1 day. This would require an exit apoapsis altitude of 71,200 km for Earth aerocapture or 
an apoapsis altitude of 33,800 km for Mars aerocapture. The data do not provide any indication of 
what the guidance scheme would be, but they do indicate a bound for the amount of control that is 
required. If the target vacuum periapsis is in the range defined by the -15° and 15° angle of attack 
curves, then there is an angle of attack within that range which would yield any desired exit apoapsis 
altitude. A successful guidance scheme would have to vary the angle of attack to continually adjust 
the trajectory during the aeropass. 

A constraint on the allowable range of target vacuum periapsis altitudes and angles of attack 
is the acceleration load experienced by the vehicle. The points at which the acceleration loads exceed 
five Earth "g's" (five times the acceleration due to gravity at the Earth's surface) are shown on the 
graphs where applicable. These g-limits decrease the range of target periapsis altitudes that can be 
used. These limits were plotted assuming constant angle of attack trajectories. In section III, it was 
shown that the peak acceleration load occurs just before the vehicle reaches periapsis. Therefore, it 
is possible to avoid exceeding the acceleration load limit by using a variable angle of attack trajec- 
tory which allows the use of a higher target periapsis altitude for the atmospheric entry. This was 
basically the approach used by Vinh and Lu that was described in section III. 

Another factor that affects the range of allowable target vacuum periapsis altitudes is the 
vehicle ballistic coefficient. The data in figures 4-5 through 4-11 show that for a given angle of 
attack, it is necessary to decrease the target periapsis altitude as ballistic coefficient increases. An 
increase in the ballistic coefficient also results in a shift of the acceleration load limit boundary. This 
means that for a given target periapsis altitude, the peak acceleration load decreases as ballistic 
coefficient increases. Unfortunately, other vehicle loads increase as ballistic coefficient increases, so 
a tradeoff must be made in selecting a vehicle ballistic coefficient. The effect of ballistic coefficient on 
the vehicle loads was investigated and results are shown in the next set of figures. 


2. Vehicle Loads 


Figures 4-12 through 4-18 show the vehicle loads at various angles of attack and ballistic 
coefficients. As was the case for the previous set of data, these data were generated assuming 
constant angle of attack trajectories. The loads for variable angle of attack trajectories would usually 
be less. The vehicle loads that were investigated included the maximum acceleration load, £ max ; the 
maximum stagnation point temperature, 7^^; and the maximum dynamic pressure, <? max . Where it 

is necessary, the 5-g acceleration limit and the maximum temperature limit of current reusable 
thermal protection systems (1,645 K) are indicated on the graphs. No limit for the maximum dynamic 
pressure is indicated since that is usually a function of vehicle design. 

The loads in figures 4-12 and 4-13 are based on an assumed exit apoapsis altitude range 
from 120 to 1,000 km. In reality, the loads decrease as exit apoapsis altitude increases, but for the 
range of exit apoapsis altitudes considered, the decrease is usually not significant. The data for the 
aerocapture applications, figures 4-14 through 4-18, have a slightly different format. A range of 
vehicle loads are depicted for a range of exit apoapsis altitudes, from the top of the atmosphere to 
about 100,000 km. The loads are shown for three angles of attack and three ballistic coefficients. The 
maximum load indicated at the top of each bar corresponds to a low exit apoapsis altitude, and the 
maximum load indicated by the bottom of each bar corresponds to a high exit apoapsis altitude. 
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Figure 4-13. Lunar return aerobraking vehicle loads. 
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Figure 4-15. Earth aerocapture vehicle loads, C 3 = 50 km 2 /s 2 . 
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Figure 4-16. Mars aerocapture vehicle loads, C 3 = 10 km 2 /s 2 . 
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The data in figures 4-12 through 4-18 show the importance of careful sizing of an aerobrake. 
The ballistic coefficient of the vehicle has a strong influence on the loads. As the ballistic coefficient 
increases (i.e., aerobrake diameter decreases), the maximum acceleration load decreases slightly, 
but the maximum stagnation point temperature and maximum dynamic pressure can increase 
significantly. The angle of attack can also have an impact on the vehicle loads. As the angle of attack 
becomes more positive (increasing lift downward) the maximum loads decrease. The sometimes 
counteracting effects of increasing ballistic coefficient and angle of attack on the vehicle loads show 
why it is necessary to consider both the vehicle design and the trajectory design in the process of 
sizing an aerobrake. 

The two previous sets of data which were generated using constant angle of attack trajec- 
tories indicate that to minimize the vehicle loads it would be desirable to fly with lift directed down- 
ward at the highest possible L/D. This would result in the nearly grazing trajectory which was de- 
fined earlier [1], Unfortunately, even at low L/D's, at positive angles of attack the exit apoapsis is 
too sensitive to the target vacuum periapsis. This is indicated by the almost vertical curves on the 
graphs in figures 4-5 through 4-11. One way to decrease the sensitivity to vacuum periapsis altitude 
and also reduce vehicle loads is to use a variable angle of attack trajectory with a lift profile similar 
to the one developed by Vinh and Lu [22]. 


3. Comparison of Constant an d Variable Angle of Attack Trajectories 

Comparisons of the trajectories and maximum vehicle loads for constant and variable angle 
of attack trajectories are shown in figures 4-19 through 4-25. Due to the limitations of the computer 
programs that were available, it was not possible to find the optimum lift profile for each aerobrake 
mission, therefore, a lift profile that approximated the one developed by Vinh and Lu [22], and de- 
scribed in the previous section, was used. The variable lift profile used a simple two-step approxi- 
mation to Vinh and Lu's solution. Maximum upward lift (negative angle of attack) was used for the 
descending leg of the aeropass, then maximum downward lift was used for the ascending leg of the 
aeropass. The switch from negative to positive angle of attack was performed at the minimum 
altitude of the trajectory and was assumed to be almost instantaneous. 

Each of the figures in this set of data show the altitude profile during the aeropass as a func- 
tion of time for three constant angle of attack trajectories and the variable angle of attack trajectory. 
The data listed with each figure show several trajectory parameters and the peak vehicle loads for 
each trajectory. The variable angle of attack trajectory is designated as "-/+15.0," which signifies an 

angle of attack of —15° during the descending leg of the aeropass and +15° during the ascending leg 
of the aeropass. 

The trajectories are typical for each application. The exit apoapsis altitude for the GEO return 
and lunar return trajectories is representative of a return to an LEO space station orbit altitude. The 
exit apoapsis altitude for the Earth aerocapture missions represents the apoapsis altitude of an 
elliptical orbit with a period of 24 hours and a periapsis altitude of about 500 km. The exit apoapsis 
altitude of the Mars aerocapture missions represents the apoapsis altitude of an elliptical orbit with 
a period of one Mars solar day and a periapsis altitude of 250 km. The use of a high exit apoapsis for 
aerocapture missions reduces the peak vehicle loads. 

A variable angle of attack lift profile provides a compromise between the disadvantages of 
constant upward and constant downward lift profiles. Trajectories flown with constant lift upward 
(negative angle of attack) are less sensitive to target vacuum periapsis altitude but can have high 
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ALTITUDE (KM) 


GEO RETURN AEROBRAKING TRAJECTORIES 

BALLISTIC COEF. = 100 KG/M“2 EXIT APOGEE ALT. = 500 KM 



TIME (SEC.) 


ALPHA 

TARGET 

PERIAPSIS 

EXIT 

AP0APSIS 

ENTRY 

VELOCITY 

(DEG) 

(KM) 

(KM) 

(M/S) 

-15.0 

59.9800 

499.95 

9899.85 

0.0 

74.1400 

499.90 

9899.68 

15.0 

80.5563 

500.57 

9899.60 

-/ + 15.0 

73.9357 

499 . 80 

9899.68 


MIN 

MAX 

MAX 

MAX 

ALT 

G-LOAD 

HEAT 

DYN 



RATE 

PRESS 

(KM) 

(9) 

( W/CM* *2) 

(N/M**2) 

69.21 

3.91 

37.20 

3735.91 

73.25 

2.28 

29.55 

2231.11 

77.01 

1.40 

23.78 

1338.67 

76.45 

1.58 

26.13 

1510.07 


Figure 4-19. GEO return aerobrake trajectories. 
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ALTITUDE (KM) 


LUNAR RETURN AEROBRAKING TRAJECTORIES 



ALPHA 

TARGET 

EXIT 

ENTRY 

MIN 

MAX 

MAX 

MAX 


PERIAPSIS 

APOAPSIS 

VELOCITY 

ALT 

G-LOAD 

HEAT 

DYN 

(DEG) 






RATE 

PRESS 

(KM) 

(KM) 

(M/S) 

(KM) 

(9) 

(W/CM* *2) 

(N/M* *2) 

-15.0 

54.5820 

499.93 

10574.89 

66.79 

5.99 

51.68 

5728.59 

0.0 

72.1521 

500.08 

10574.32 

71.32 

3.23 

39.82 

3170.00 

15.0 

79.3996 

502.46 

10574.08 

75.59 

1.92 

31.48 

1832.77 

-/+ 15.0 

68.5790 

500.49 

10574.43 

73.05 

2.89 

39.31 

2764.95 


Figure 4-20. Lunar return aerobrake trajectories. 
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ALTITUDE (KM) 


EARTH AEROCAPTURE TRAJECTORIES 



ALPHA 

TARGET 

PERIAPSIS 

EXIT 

APOAPSIS 

ENTRY 

VELOCITY 

(DEG) 

(KM) 

(KM) 

(M/S) 

-15.0 

75.8992 

71225.75 

11104.50 

0.0 

78 . 7908 

71227.04 

11104.40 

15.0 

81.1538 

71226.40 

11104.31 

-/+ 15.0 

77.0565 

71225.33 

11104.46 


MIN 

MAX 

MAX 

MAX 

ALT 

G-LOAD 

HEAT 

DYN 



RATE 

PRESS 

(KM) 

(9) 

(W/CM**2) 

(N/M**2) 

77.87 

1.70 

35.53 

1625.45 

78.60 

1.46 

32.98 

1432.51 

79.47 

1.30 

31.20 

1240.04 

78.74 

1.49 

33.61 

1428.68 


Figure 4-21. Earth aerocapture trajectories, C 3 = 10 km 2 /s 2 . 




ALTITUDE (KM) 


EARTH AEROCAPTURE TRAJECTORIES 
BALLISTIC COEF. = 100 KG/M**2 

C3 = 50 KM**2 / S**2 EXIT APOAPSIS ALT. = 71225 KM 



ALPHA 

TARGET 

EXIT 

ENTRY 

MIN 

MAX 

MAX 

MAX 


PERIAPSIS 

APOAPSIS 

VELOCITY 

ALT 

G-LOAD 

HEAT 

DYN 

(DEG) 






RATE 

PRESS 

(KM) 

(KM) 

(M/S) 

(KM) 

(9) 

(W/CM**2) 

(N/M**2) 

-15.0 

64.0276 

71225.80 

12727.35 

70.06 

6.02 

82.38 

5752.05 

0.0 

72.4564 

71226.57 

12726.99 

72.12 

4.38 

71.42 

4290.82 

15.0 

77.7764 

71214.72 

12726.76 

74.43 

3.33 

61.57 

3179.57 

-/+ 15.0 

68.5179 

71224.53 

12727.16 

72.55 

4.41 

73.45 

4215.65 


Figure 4-22. Earth aerocapture trajectories, C3 = 50 km 2 /s 2 . 
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ALTITUDE (KM) 


MARS AEROCAPTURE TRAJECTORIES 
BALLISTIC COEF. = 100 KG/M“2 

C3= 10 KM**2 / S**2 EXIT APOAPSIS ALT. = 33793 KM 



TIME (SEC.) 


ALPHA 

TARGET 

PERIAPSIS 

EXIT 

APOAPSIS 

ENTRY 

VELOCITY 

(DEG) 

(KM) 

(KM) 

(M/S) 

-15.0 

45,0521 

33793.52 

6073.00 

0.0 

51.2504 

33793.43 

6073.22 

15.0 

55.8255 

33791.74 

6073.39 

-/+ 15.0 

48.3351 

33793.65 

6073.12 


MIN 

MAX 

MAX 


MAX 

ALT 

G-LOAD 

HEAT 

RATE 


DYN 

PRESS 

(KM) 

(9) 

(W/CM 

1 ** 2 ) 

(N/M**2) 

49.13 

1.46 


8.30 

1397.08 

50.78 

1.15 


7.41 

1128.48 

52.60 

0.94 


6.73 

893.88 

51.36 

1.14 


7.45 

1093.53 


Figure 4-23. Mars aerocapture trajectories, C 3 = 10 km 2 /s 2 . 
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MARS AEROCAPTURE TRAJECTORIES 
BALLISTIC COEF.= 100 KG/M“2 
C3= 25 KM"2 / S**2 


EXIT APOAPSIS ALT. = 33793 KM 



ALPHA 

TARGET 

EXIT 

ENTRY 


PERIAPSIS 

APOAPSIS 

VELOCITY 

(DEG) 

(KM) 

(KM) 

(M/S) 

-15.0 

34.2682 

33792.97 

7234.13 

0.0 

46 . 0666 

33793.06 

7234.63 

15.0 

53.1162 

33788.04 

7234.92 

-/+ 15.0 

42.3419 

33795.40 

7234.47 


MIN 

MAX 

MAX 

MAX 

ALT 

G-LOAD 

HEAT 

DYN 



RATE 

PRESS 

(KM) 

(9) 

(W/CM**2) 

(N/M**2) 

42.56 

4.03 

18.30 

3850.82 

45.54 

2.76 

15.44 

2706.06 

48.63 

1.96 

13.30 

1869.15 

46.98 

2.55 

15.55 

2436.99 


Figure 4-24. Mars aerocapture trajectories, C 3 = 25 km 2 /s 2 . 
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ALTITUDE (KM) 


MARS AEROCAPTURE TRAJECTORIES 
BALUSTIC COEF.= 100 KG/M**2 
C3= 40 KM**2 / S“2 


EXIT APOAPSIS ALT. = 33793 KM 



TIME (SEC.) 


ALPHA 

TARGET 

PERIAPSIS 

EXIT 

APOAPSIS 

ENTRY 

VELOCITY 

(DEG) 

(KM) 

(KM) 

(H/S) 

-15.0 

27.0874 

33793.22 

8229 . 39 

0.0 

43.4275 

33791.81 

8230.13 

15.0 

51.9733 

33810.53 

8230.51 

-/+ 15.0 

36.6991 

33793.68 

8229.82 


MIN 

MAX 

MAX 

MAX 

ALT 

G-L0AD 

HEAT 

DYN 



RATE 

PRESS 

(KM) 

(9) 

( W/CM* *2) 

(N/M**2) 

39.15 

6.93 

30.52 

6628.69 

42.96 

4.41 

25.42 

4320.50 

46.90 

2.93 

21.08 

2800.86 

43.74 

4.71 

26.74 

4500.52 


Figure 4-25. Mars aerocapture trajectories, C 3 = 40 km 2 /s 2 . 
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peak loads. Trajectories flown with constant lift downward (positive angle of attack) have lower 
peak loads but are very sensitive to target vacuum periapsis altitude. The variable angle of attack 
profile results in a minimum altitude between that of a constant lift up or constant down profile. This 
results in peak loads which are between the two extremes but generally closer to the constant lift 
down trajectory and a sensitivity to target vacuum periapsis that is similar to a zero lift trajectory 
(constant angle of attack equal to zero). 


All ot the trajectories in figures 4-19 through 4-25 were generated using the standard 
atmosphere models described earlier. Since the density profiles were smooth functions of altitude, 
the lift profiles, whether constant or variable, represent "nominal" or baseline cases. In reality, the 
lift would have to vary from the nominal to correct for guidance errors and atmospheric variations. 
Varying the angle of attack allows the trajectory to be adjusted during the aeropass. The benefits of 
a variable angle of attack can be seen by looking at the effect of density variations on aerobrake 
trajectories and the range of allowable target vacuum periapsis altitudes. 


4 - Definition of Target Vacuum Periansis Tolerance 


As stated earlier, the actual guidance scheme that would allow a vehicle to achieve an exact 
exit condition was not investigated. Rather, given an assumed vehicle ballistic coefficient and angle 
of attack range, the limits within which the guidance scheme would operate were defined. Figures 
4-26 through 4-28 show the range of allowable target vacuum periapsis altitudes for an aerobraked 
vehicle at various atmospheric densities for each aerobrake application. The vehicle was assumed to 
have a ballistic coefficient of 100 kg/m 2 and an angle of attack control range from -15° to +15°. Tra- 
jectories were analyzed using atmosphere models corresponding to the mean density profile, 20 per- 
cent above the mean and 20 percent below the mean. In reality the deviation of the density from the 

mean is not constant with altitude so this data provides conservative boundaries of the density 
variation. J 

Each bar in these figures represents a range of target vacuum periapsis altitudes which allow 
an exit apoapsis within the designated acceptable range without exceeding a maximum acceleration 
load of five Earth g's. For GEO and lunar return (fig. 4-26), the designated range of exit apoapsis 
altitudes is from 130 to 1,000 km. For the aerocapture applications (figs. 4-27 and 4-28) the ranee 
is from 500 to 100,000 km. 6 


The range of altitude for which the three bars overlap define the target vacuum periapsis 
altitudes for which the vehicle could achieve a desired exit apoapsis as long as the density is within 
20 percent of the mean. This range of altitudes will be defined as the target vacuum periapsis altitude 
tolerance. The size of this tolerance is shown for each aerobrake application. It defines the accuracy 
to which the vehicle must be targeted prior to atmospheric entry. Once the vehicle enters the 
atmosphere, the trajectory would be adjusted by varying the angle of attack within the control range. 

To find the limits of the target vacuum periapsis altitude without defining a guidance scheme, 
a very conservative approach was used for adjusting the trajectory in the atmosphere. Only one 
angle of attack adjustment was used. The entry angle of attack was held constant until the vehicle 
reached minimum altitude, then the angle of attack was adjusted to either maximum lift down or 
maximum lift up for the ascending leg of the trajectory. This single large adjustment reveals the 
limits of what a real guidance scheme could achieve with many small adjustments. 
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VACUUM PERIGEE ALT. (KM) 


EXIT APOGEE ALT. RANGE = 1 30 - 1 ,000 KM 
KEY: 0 -20% DENSITY VARIATION 
□ MEAN DENSITY PROFILE 



ENTRY ALPHA = -15 DEG ENTRY ALPHA = -15 DEG 

EXIT ALPHA = +15 DEG EXIT ALPHA = +15 DEG 


Figure 4-26. Effect of density variations on target vacuum periapsis tolerance for 
* GEO and lunar return aerobraking. 


EXIT APOGEE ALT. RANGE = 500 - 100,000 KM 
KEY: 0 -20% DENSITY VARIATION 
□ MEAN DENSITY PROFILE 


£13 +20% DENSITY VARIATION 



Figure 4-27. Effect of density variations on target vacuum periapsis tolerance 
* for Earth aerocapture. 
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EXIT APOAPSIS ALT. RANGE = 500 - 100,000 KM 
KEY: £3 -20% DENSITY VARIATION 
□ MEAN DENSITY PROFILE 



ENTRY ALPHA =-15 DEG ENTRY ALPHA =-15 DEG ENTRY ALPHA = 0 DEG 

EXIT ALPHA = +15 DEG EXIT ALPHA = +15 DEG EXIT ALPHA = +15 DEG 

Figure 4-28. Effect of density variations on target vacuum periapsis tolerance 

for Mars aerocapture. 


In agreement with the results of Vinh and Lu [22], it was found that the vehicle loads were 
reduced if maximum upward lift (angle of attack equal to -15°) was used during the descending leg of 
the trajectory. In most cases, this entry angle of attack also resulted in the largest target vacuum 
periapsis altitude tolerance. The exceptions were the aerocapture applications at high C 3 's where it 
was found that a ballistic entry (angle of attack equal to 0°) yielded the largest target vacuum 
periapsis tolerance. 


V. SUMMARY 


This report has investigated the application of aerobraking for orbital transfers from geosyn- 
chronous and lunar orbit to low Earth orbit and aerocapture at Earth and Mars. The objective was to 
show that low L/D aerobrakes which use angle of attack variation for control could be used for each 
of these applications. The investigation began with a study of the design considerations that are 
necessary in the process of designing and sizing an aerobrake for a particular application. Since the 
aerobrake vehicle design process is interrelated to the trajectory design, a review of several formu- 
lations for obtaining optimum aerobrake trajectories was presented. Finally, the relationship 
between aerobrake design factors and trajectory characteristics was investigated further using 
numerical simulations to analyze aerobrake trajectories for each application. 
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The aerodynamic characteristics of aerobrakes are in many ways unique from most other 
atmospheric vehicles. The one characteristic that sets aerobrakes apart is that they are designed 
primarily to utilize the drag force rather than the lift force. Lift is used only as a control for changing 
the lateral direction of flight or the rate of altitude change. In many aerobrake applications, very little 
lift is needed. Each of the aerobrake applications that were investigated could be implemented with a 
maximum L/D of only 0.3 using a spherical aerobrake. It was shown that the L/D is closelyrelated to 
the ballistic coefficient in determining the atmospheric trajectory of an aerobraked vehicle. The proper 
combination of L/D and ballistic coefficient allows the vehicle to achieve sufficient depth of penetra- 
tion into the atmosphere for aerodynamic control. This combination is determined by the vehicle 
acceleration and heating load limits. 


Due to the unique aerodynamic characteristics of aerobraked vehicles, special attention must 
be paid to the aerodynamic stability of the vehicle. The conditions for static longitudinal stability 
were investigated. It was shown that the location of the vehicle center of mass must be in front ot 
the neutral point so that C ma is less than zero. The neutral point was defined as the point on the 
vehicle trim line where C ma equals zero. The location of the neutral point depends on the aero- 
dynamic coefficients and the diameter of the aerobrake. For a spherical aerobrake, the distance from 
the front of the aerobrake to the neutral point increases as the angle of attack decreases or the 
aerobrake diameter increases. 


Several other aerobrake design factors were also described. One of the most important 
design considerations is the flow behind the aerobrake. The aerobrake must be sized so that flow 
does not impinge on the vehicle behind the aerobrake. For angles of attack between 5 and 15 , the 
flow impingement angle ranges from 20° to 35° of inward deflection behind the aerobrake. The flow 
impingement places limits on the size of the aerobrake payload since the aerobrake diameter must 
often be balanced against the capabilities of the vehicle to withstand thermal and aerodynamic loads. 
As indicated previously, the effect of aerobrake diameter on the vehicle design must be carefully 
considered since the ballistic coefficient of the vehicle has a large influence on the atmospheric tra- 
jectory, which in turn determines the magnitudes of the loads. 


One of the most critical vehicle loads that influences the design of an aerobrake is the aero- 
dynamic heating. The heating loads determine to a large extent the size of the aerobrake. The mag- 
nitude of the heating loads that are allowed is determined primarily by the capabilities of the TPS. 
Several types of TPS's were reviewed, but due to the large heat loads that can be experienced, it 
appears that reusable surface insulation is the best choice. The aerodynamic heating environment for 
an aerobrake is not yet clearly understood. There is considerable uncertainty concerning the 
existence and effects of nonequilibrium aerodynamic heating during an aerobraking maneuver. One ot 
the primary objectives of the planned AFE is to investigate this phenomenon. 


The optimization of aerobrake trajectories depends largely on how optimum is defined. There 
are many performance indices that can be used as criteria for defining an optimum trajectory. In many 
cases the optimum trajectories, as measured by different performance indices, are very different. It 
was found that performance indices which result in "nearly grazing" optimal trajectories provided the 
best compromise between the often conflicting goals of minimizing the vehicle propulsive require- 
ments and minimizing the aerodynamic healing. To implement an optimum aerobraking trajectory, it 
is necessary to select a vehicle control strategy. The performance indices and allowable trajectory 
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controls are closely related to or constrained by the vehicle design. This indicates that there is an 
interrelationship between vehicle design and trajectory design. Roll control, pitch control (angle of 
attack variation), and drag control were discussed. 

The investigation ot the relationship between vehicle and trajectory design and the potential 
for using angle of attack variation for control was continued using numerical simulations to analyze 
trajectories for each aerobrake application. The data that were presented showed the sensitivity of 
the trajectory to vehicle angle of attack, ballistic coefficient, and entry conditions. A conservative 
approach was used in modeling the atmospheric density variations and aerobrake angle of attack 
control profile. This provided a representative boundary of the effects of atmospheric density 
variations on the trajectories and the capabilities of low L/D aerobrakes across a wide range of 
trajectory conditions. c 


The data that were developed showed that aerobrake trajectories are very sensitive to the 
atmospheric entry conditions and the density variations within the atmosphere. This implies that a 
very responsive control scheme is required. One potential advantage in the responsiveness of pitch 
control over roll control is that there is no lateral component of the lift force to be balanced in order to 
maintain a particular orbital plane. Data were also developed to illustrate the relationship between 
vehicle design and trajectory design. The sensitivity of several vehicle loads to angle of attack and 
ballistic coefficient was investigated. The data showed that in some applications the acceleration and 
heating constraints may limit the range of allowable aerodynamic controls. To reduce the peak loads 
it may be necessary to constrain the magnitude and direction of the lift force to limit the depth of 
penetration into the atmosphere. 

The feasibility of using low lift to drag aerobrakes for the selected applications was investi- 
gated further by determining the range of target vacuum periapsis altitudes for which an aerobraked 
vehicle could achieve desired atmospheric conditions. A simple two-phase angle of attack control 
profile was used to provide an approximation of the optimum lift profile to reduce vehicle loads and 
provide an approximation of the control limits for low 1VD aerobrakes. The target vacuum periapsis 
tolerances that were defined indicate that low L/D aerobrakes have sufficient control capability to 
perform each of the aerobrake applications. 

The aerobrake trajectory analysis that was performed for this report provides only a shallow 
penetration (aerobrake pun) into some of many problems that will require thorough investigation if 
aerobraking is to become a useful operation for space transportation vehicles. Nevertheless, the 
results in this report have shown that a wide range of aerobrake applications are feasible for vehicles 
with low L/D aerobrakes using angle of attack variation for control. 
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